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We study the electric current in the non-equilibrium Kondo model at zero magnetic field, using 
real-time perturbation theory in the Schwinger-Keldysh formulation. We show that the perturbative 
coefficients to all orders have a finite limit at large switch-on time (to — > — oo), and we give a 
prescription for general operators to give finite coefficients in this limit. We explain how this is 
related to the fact that the leads play the role of thermal baths and allow relaxation to occur and 
the steady state to form. This proves perturbatively that a steady state is reached in the Schwinger- 
Keldysh formulation, and specifies which operators correspond to quantities that have a well-defined 
value in the steady state. Then, we show that the steady state can be described by a special type of 
density matrix (related to Hershfield's conjecture for the particular example of the non-equilibrium 
Kondo model.) In the second part of the paper we perform a renormalization-group analysis of 
the perturbative series. We give a general argument that strongly suggests that the perturbative 
series of any average in the steady state satisfies the equilibrium Callan-Symanzik equations, and 
show in detail how it works to one-loop order for the electric current operator inside any average. 
We finally compute to two loops order the average of the electric current in the steady state, and 
perform a renormalization-group improvement. From this, we give a universal prescription, valid in 
the perturbative regime, for comparing the effect of the electric current to that of the temperature 
on the "Kondo cloud". 



INTRODUCTION AND DISCUSSION 



The description of an out-of-equilibrium strongly correlated system is a long standing problem. Even in the 
simplest case where the system is in a steady state and its properties no longer change with time, the usual 
formalism of quantum statistical mechanics is inadequate. Theoretical understanding of such systems became all 
the more pressing with the recent spectacular progress in nanotechnology, which has made it possible to study 
the Kondo impurity, one of the best understood strongly correlated systems, in out-of-equilibrium conditions. 

The Kondo impurity was realized experimentally as a quantum dot, a tiny island of electron liquid, attached via 
two tunnel junctions to leads (baths or reservoirs of electrons) held at different electric (or chemical) potentials. 
This set-up allows an electric current to flow across the dot, and measurements of the current were carried out 
as a function of the potential difference V, the temperature T and the magnetic field B 

When the dot carries a net spin in the Coulomb blockade regime, it can be modeled by a Kondo Hamiltonian 
with two channels a = 1,2, corresponding to the two leads, to which the spin of the dot, S, couples Q. The 
resonant tunneling through the dot (elastic co-tunneling) allows the electrons from each bath to jump on the 
dot and back to the same bath, leading to the formation of Kondo resonance around the Fermi level fi a in 
each lead. Further, electrons from one bath can jump on the dot and onto the other bath, giving "off-diagonal 
coupling" of the two channels to each other. With the matrix of couplings J a , a ' and at zero magnetic field, the 
Hamiltonian is 



a k,a a > a ' k,k',a,a' 

Here a denotes the spin index a = ±1/2, and S is in the spin-1/2 representation. 



The process corresponding to off-diagonal coupling induces a current when the baths are held at nonzero 
potential difference, V = [ii — The development of the Kondo resonance as temperature is lowered enables 
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FIG. 1: (Color online) (A) The system. (B) The system unfolded. 



the system to overcome the Coulomb blockade, producing a significant increase of the conductance. The unitarity 
limit is reached as T — > 0, V — > 0. 

As we are interested in the universal properties of the system, we shall consider the model in the range T, V <C 
D a , where D a = D are the bandwidths of the leads, each lead being considered a very large conductor (the 
bandwidths can be assumed to be the same for both channels). We are allowed therefore to carry out the 
standard steps (linearizing around the Fermi level, keeping only the s-wave component in the expansion of c Q g 
in spherical modes), to obtain a representation of each lead as a free electron gas on the half line consisting 
of left and right movers i/j aj L(x),ip a! ii(x),x < 0, interacting with the impurity localized at x = 0. It will be 
convenient for us to "unfold" the baths, making left and right movers on the half line into right movers on the 
full line (defining ip a ,B,(%) = ^oi,l{—x),x > 0). See Figure 1. 

The field-theoretic Hamiltonian is then: 

/°° V f°° 

dxil>l{x)di) a {x) + - dx(V#2-w)+ J2 J <*,<x'i>l(0)^<x'(0)-S 

= H + VH Z +H! (2) 
where we work in units such that vf = 1. We also denote, 

1 f°° 

H z = - dx (i>\i>2 - i\)\i\>\) 



2 . 
and 

Hl= Jc t ,a>1>U0W«(0)-& 



The coupling of the baths to the quantum dot is parametrized by the Hermitian matrix 

J a ,af = * ( ^ ^ ) (3) 
\ / a, a' 

where and A are real (the factor tt is introduced for later convenience). It is possible to diagonalize the matrix 
of coupling J QjCt ' by a change of basis in the channel space. In the simple situation of a single level quantum 
dot, described by the Anderson model, coupled to identical leads one naturally finds the relation A^ = A. This 
makes the matrix of coupling constants degenerate, one of the eigenvalues being zero. After diagonalizing it, 
the terms representing the interaction with S become that of a decoupled free fermion and a one-channel Kondo 
model. More generally, for A^ =/= A, diagonalizing the matrix of couplings gives the interaction term of the usual 
2-channel Kondo model. Out of equilibrium, for V ^ 0, diagonalizing the matrix of couplings does not lead to 
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a simplification of the problem, since the out-of-equilibrium term VH Z induces extra coupling between the new 
fermion fields: it is not invariant under change of basis. We will not perform this diagonalization here, in order 
to keep the term VH Z simple. Also, we will consider the general case A^ 7^ A. 

In this paper, we are interested in studying the electric current as function of voltage and temperature in the 
steady state of this model. As in the usual Schwinger-Keldysh formulation of a non-equilibrium steady state 
[IS EH , we will imagine coupling the dot to the leads at some time t in the past when the system is in a 
thermal equilibrium state, then allowing the system to evolve till t = when the current is evaluated. One 
expects that after a transitory regime, as to * —00, the system will relax into a steady state with a constant 
current flowing through the dot. 

Many questions arise concerning this formulation. Probably the most obvious one is: Is the model sufficient to 
describe the establishment of a steady state current? Or does one have to contemplate additional relaxation 
mechanisms (certainly present in actual experiments) to absorb the continuous flow of energy of the electrons 
moving from the higher Fermi-level lead to the lower one? In the framework of real-time perturbation theory, a 
related (but not equivalent) question that one can answer is whether or not the "infrared" limit to — > —00 exists 
for the integrals representing the perturbative coefficients in an expansion in A and A,;. We develop the real-time 
Keldysh perturbation theory (in some ways similar to [H and [l7j])> and use it to establish the convergence of 
every term of the perturbative series as the switch-on time to is sent to minus infinity. 

This result is highly non-trivial. Real-time perturbation theory often gives divergences as the switch-on time is 
sent to minus infinity, unless a good relaxation phenomenon is included in the model. In equilibrium, this issue 
can be easily overcome: an infrared-divergent real-time perturbation theory for a system in equilibrium only 
means that the particular model we are considering does not have the proper relaxation mechanism. But with an 
additional interaction, however small, representing a proper relaxation mechanism, the system will have infrared 
convergent real-time perturbation theory. Then, for equilibrium models, this can be equivalently described by 
the always-infrared-convergent "imaginary time" perturbation theory (where the integrals in imaginary time 
are on a finite interval) coming from the description of the model using its equilibrium density matrix. There, 
the additional interaction representing relaxation can be sent to zero from the beginning. 

Out of equilibrium, however, there is no a priori steady-state density matrix description of steady quantities. 
Thus, if the real-time perturbation theory is infrared divergent, there is no simple way to describe steady state 
physics. This is very natural: in contrast to the equilibrium case, the relaxation mechanism is used not only 
to reach the steady state, but also to form it, since we need a continuous absorption of energy. Hence, such 
infrared divergences are far more pathological. If the limit of large negative switch-on time exists order by order, 
this superficially seems a good indication that the steady state is reached and that the model indeed describes 
the steady state (although, strictly speaking, one would still have to analyse the possible non-perturbative 
contributions if the perturbative series is asymptotic). But if no such limit exists, certainly more questions 
arise: does the model reach a steady state non-perturbatively (that is, divergences are an artifact of perturbation 
theory), or are other interactions necessary? These questions of course appeared in the literature before (see 
for instance [H, EH]). An interesting example is the case studied in [12], where it was shown that when the 
model we are discussing is put into a magnetic field, the limit to ~~ * — °o and A, Ad — > do not commute, 
leading to divergencies in the perturbation theory as to — > —00. Assuming that a steady state exists in the 
model, the correct result as to — > —00 should then be non-perturbative, and it was partially evaluated under 
this assumption (the "zeroth order" was evaluated), without external thermal reservoir coupled to the dot. 

But the convergence of the perturbative series or the technical "way around" its divergencies as just described 
above do not guarantee that the results are describing the correct physical steady state. Indeed, if the model 
is believed to have a steady state, perturbatively or not, two questions should still be answered: Is there an 
element in the model playing the role of a good thermal reservoir to sustain the correct steady state in the 
model? And if not, is there a guarantee that a coupling to an external reservoir would not have an important 
effect? 

These last two questions have more bearing than it may seem. In the usual Schwinger-Keldysh formulation, 
one does not assume any exchange with a thermal bath while the system is evolved: one starts with a thermal 
equilibrium state, then turns on the coupling to the dot and lets the system evolve without thermal bath. This 
is certainly not the true physical situation; in fact, understanding how a thermal bath affects the evolution 
of a quantum system was amongst the main points of the study of Caldeira and Leggett We carry out 
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the real-time perturbation theory both in equilibrium (that is, in the case V = 0; note that the real-time 
formalism includes a transitory non-equilibrium region where the system relaxes to equilibrium) and out of 
equilibrium (V ^ 0), and show that it is the same phenomenon that makes the perturbative series infrared 
convergent in equilibrium and out of equilibrium. This phenomenon is a factorization at large time separation 
of the correlation functions involved in the perturbative coefficients; interpreting the integrals over time defining 
these coefficients in terms of physical processes, this signals a decoherence in time induced by the leads and 
suggests that the leads are good thermal baths. The convergent expression in equilibrium is indeed the right 
equilibrium density matrix, confirming that the leads themselves play the role of thermal baths for the dot 
degrees of freedom. The convergent expression out of equilibrium should then be "the right" steady state. Note 
that the factorization signaling decoherence in time occurs because the Hamiltonian for the leads is conformally 
invariant: indeed, in conformal field theory, a large time separation is a large distance separation, which, by 
locality, gives rise to factorization. Physically, this occurs because the separation between the energy levels of 
the baths is much smaller than all other scales in the problem and essentially energy-independent (although, 
as will become clear from our investigation, these two conditions may not be sufficient). Note that these ideas 
are not entirely new: in Q the results of Caldeira and Leggett for constructing a thermal bath coupled to a 
quantum mechanical system were re-interpreted as coupling a bulk conformal field theory in a disk to degrees 
of freedom on the boundary. Our results generalize this to the behavior of the impurity in the Kondo model. 
We are currently investigating how this can be further generalized (for example, what the general properties of 
the conformal field theory should be). 

Our proof also allows us to describe the steady-state physics in terms of a "steady state density matrix," as 
conjectured by Hershfield 0]. The essential difference between the usual density matrix and the steady-state 
density matrix can be seen as a non-locality in the latter which captures the build-up of the steady state. 

Note that a proof of convergence to all order was developed in 3] for a non-equilibrium free boson model 
with boundary interactions, but the arguments there were very model-dependent (and do not apply to the 
non-equilibrium Kondo model) and quite different from ours. In particular, our arguments have a much deeper 
physical meaning and scope, as explained above. 

Let us stress here that it is quite important to know that no external bath is required for reaching a non- 
equilibrium steady state in an impurity model. Indeed, this means that we can use the real-time formalism 
without addition of a coupling to an external bath in order to 1) study the perturbative scaling properties of the 
model as we did in this paper for the Kondo model (see below), 2) construct more or less explicitly, without 
strong assumption, the steady state as an eigenstate of the Hamiltonian, for instance using the the formalism 
of the "steady-state density matrix" ; this eventually can give access to the infrared behavior of the model and 
to the integrability properties of the steady state (work in progress). The real-time formalism might not be 
the easiest way of trying to obtain this understanding, but it is probably the clearest, as it is the most closely 
related to the actual experimental situation. 

Other questions that need to be addressed in the study of this out-of-equilibrium steady state are: How will the 
Kondo effect, the quenching of the impurity spin as the temperature is lowered below the Kondo scale Tk, evolve 
in the presence of a current? Will new scales make their appearance? Which quantities are universal? To what 
extent the powerful ideas of the Renormalization Group (RG) apply there? Many interesting attempts were 
carried out, mainly perturbatively, to understand the flow of couplings as the cut-off (band width) D is reduced 
0, H, H, EH EH H3| • We shall pursue a different track and study a question related to the universality features 
of the model, namely: does a limit D — > oo exist? In this limit all results are universal. We shall establish that 
such a limit exists by running the RG equations "backwards", referring to them in the field theoretic context, 
in the usual way, as the Callan-Symanzik equations. We shall deduce an out-of-equilibrium /3-function carrying 
out the calculation directly in the steady state and will show that it is the same as the equilibrium /3-function. 
This may not be too surprising since the singularity structure of the system usually does not depend on the 
state in which they are evaluated, so that the ground state and the highly excited steady state produce the 
same singularities. The finite parts of course are different. We show that only one scale arises, the Kondo 
temperature, Tk, and the current can be written in a universal form as a function of the ratios T /Tk, V/Tk, C , 
with C an additional dimcnsionlcss parameter characterizing the asymmetry between A and Xd (in other words, 
specifying the RG trajectory). We carry out the computation of the current to two- loop order and verify these 
statements explicitly. We then use the RG arguments to re-sum the leading logarithms. Our results are valid 
in the regime where both the bias voltage and the temperature are smaller than the band width, and where 
the bias voltage or the temperature is larger than the Kondo scale: V, T <C D, and Tk <C V or Tk <C T. In 
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particular, we verify that there are no divergencies at T/V — > in the perturbative results. This means that to 
two-loop order, the voltage plays the role of a good infrared cutoff. From the RG analysis, we give a universal 
prescription, valid in this regime, for comparing the effect of the electric current to the effect of the temperature 
on the destruction of the Kondo cloud. 

We also examined the effect of a local magnetic field on the dot but, as expected, were unable to show that the 
perturbation series converges in this case. We will come back to a discussion of this case, in relation with the 
results of [Hj], in the last section of this paper. 



II. FORMULATION OF THE PROBLEM AND GENERAL CONSIDERATIONS 



The Schwinger-Keldysh formulation |l5l Il6l]. 

First formulation. 

We shall be interested in the electric current that passes from lead 2 to lead 1 across the quantum dot under 
the action of the potential diffence V. It can be calculated by evaluating the average of the current operator 
J with respect to a density matrix that has evolved over sufficiently long time from the initial non-interacting 
density matrix 

Po = e~ &H * 

under the action of the full evolution operator 

S {V) {t u t 2 ) = e-'C*!-*^ = e -i(ti-t2)(H +V H g +H T ) _ ^ 



The operational meaning of this formulation is the following. The non-interacting leads are initially, say at time 
to, brought to thermal and chemical equilibrium at zero potential difference exchanging energy and particles 
with a common external reservoir at fixed temperature and chemical potential. The energy levels of lead 1 and 
lead 2 are filled up to the same energy (with thermal and particle fluctuations). 

Just after time to, they are separated from the external reservoir, then a potential difference V is applied and 
the interaction is turned on. The application of the potential V just after time to, as usual, causes a raising of 
the energy levels of lead 2 with respect to those of lead 1. For the clarity of the discussion below, it is worth 
being more precise here. One should imagine both leads having a continuum of available states from the bottom 
of their bandwidths with increasing energies (the energies grow in a continuous way for infinite leads, of course, 
so one should think about densities of states). At time to, the available states of the leads are filled up to equal 
energies. Then, just after time to, when the potential is applied, one shifts the energies of all states of lead 
2 by, say, V/2 (towards higher energies), and the energies of those of lead 1 by — V/2, without changing the 
occupations of the states. Hence, the levels of lead 2 are now filled up to a higher energy than those of lead 1. 
Since the reservoir is disconnected and the interaction is turned on, there is a current. The steady state current 
is obtained after an infinite time, which we will take to be time (that is, we will take to — * -co). 

In the equation Q , the raising of the energy levels and the turning on of the interaction strength seem instan- 
taneous and simultaneous. But one can multiply both terms VH Z and Hj by a factor that smoothly increases 
from at time to = — oo to I at time 0, for instance the factor e _/it , in order to implement a simultaneous 
adiabatic increase of both the potential and the interaction strength. Sending fi — > (the adiabatic increase 
occurring far in the past) gives the steady state. This is really what is understood in this formulation. 

It is not obvious, a priori, that this formulation represents adequately the usual experimental situation, where 
the leads and the quantum dot are always connected to a common thermal reservoir (but not a reservoir of 
electrons), even while the steady state is being reached. However, it is natural to think that the leads can 
themselves play the role of thermal reservoirs. As discussed in the introduction, this is indeed the case, and will 
be made more precise below. 
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The average of an operator O in the steady state is then given by 



ln , v Tr[^>(0,t o ) e^gW(y) O] 

{ ° )ss = t^oo Tr[e-^°] ■ (5) 



The operators act in the Hilbert space for Hq (which is a tensor product of the two-channel free massless 
fermion Hilbert space and of the impurity space) obtained by imposing asymptotically vanishing conditions for 
the fermion fields (correlation functions of fermion fields vanish at infinite distance from each other and from the 
dot). To be more accurate, we could start by taking the fermion fields on a line segment of length L containing 
the dot, with some free boundary conditions; for instance, "0 Q (— L) = tjj a {L) (this corresponds to the usual free 
boundary conditions when fermions are folded back on the half line), then send L to infinity. The steady state 
would be obtained in the region 

L- 1 « |i | 1 « V,T (6) 
where the energy scale of switch-on, |io| -1 j suffices to smear out the energy level spacing L . 
Second formulation. 

Another formulation can be given. The initial idea of this second formulation is that the current can be created 
not only by a shift of the energies of the states of leads 1 and 2 (coming from the application of an electric 
potential), but also by putting additional electrons in lead 2 and taking away electrons from lead 1. In order to 
implement this, one starts again, at time to, with the uninteracting leads, both connected to a common thermal 
and particle reservoir, and in thcmal and chemical equilibrium; but now the chemical equilibrium is not at 
potential difference 0, but rather at a potential difference —V. The initial density matrix is then 

po = e-« ff »-™') . 

This potential difference shifts towards lower energies the states of lead 2 with respect to those of lead 1 by an 
amount V. But since there is equilibrium, the states of lead 1 and lead 2 are still filled up to the same energy. 
Note that then, as compared to the first formulation at time to, there are more available sates of lead 2 and less 
of lead 1 that are filled. 



Just after time to, the reservoirs are disconnected, then the potential is set to and the interaction is turned 
on. The density matrix po then evolves with the evolution operator at zero bias voltage S(t\,t2), 

S(t 1: t 2 ) = S< v=0 >(ii,i 2 ) = e ~^-^(H + Hl ) (7) 

Putting the potential to has the effect of raising the energy levels of lead 2 with respect to those of lead 
1 by an amount V, the same effect that occurs in the first formulation just after time to when the potential 
is applied. In contrast, though, this brings us to a situation where the available states of leads 1 and 2 have 
exactly the same energies as in the first formulation at time to (that is, at potential 0), but with more states 
filled in lead 2 and less in lead 1, so that the leads are filled up to unequal energies. This indeed implements 
having put additional electrons in lead 2 and extracted electrons from lead 1. With the interaction on and the 
reservoir disconnected, a current is created. Again, after an infinite time, the steady state should be reached. 
The current is then given by 

,^ y Tr[S(0,to)e-^- v ^ S(t o ,0) j] 

{J)sS = to^oc Tr [ e -P(H -VH z)] ( g ) 

(for more general operators O, see (J32J) ) . 

If the size of the bandwidth can be sent to infinity (when evaluating quantum averages of operators that give 
finite results in this limit), then the operational description above for the second formulation is equivalent to 
that of the first formulation, since then only the Fermi energies of leads 1 and 2 matter. In particular, raising 
the energy levels or filling states with electrons are exactly the same operation in this case. 

However, there is another difference between both formulations. In the second formulation, we can now think 
about putting a factor e _/it for adiabatically increasing the interaction strength, but there is no such possibility 
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for adiabatically increasing the potential (one would have to add the term — (1 — e~ mut )VH z in the evolution 
Hamiltonian). In other words, for practical calculations, this second formulation naturally implies that the 
energy levels of lead 2 are raised (or the states are filled), with respect to those of lead 1, by an amount V 
instantaneously, and that the interaction is then turned on adiabatically. This is to be contrasted with the 
first formulation, where both the potential difference and the interaction strength were understood as being 
simultaneously increased adiabatically. 

We will show below that both formulations are equivalent. 



Symmetry currents. The Hamiltonian Hq is conformally invariant and has a large algebra of symmetries 
associated with it. It is a WZW "current algebra" of the symmetry currents (not to be confused with the 
physical current and it will be convenient to carry out many of the calculations in terms of symmetry 
currents. Introduce the following operators, 

J z = tt : (iplip2 - ip\ipi) ■ 
J x = iw(i/jlaipi - ^\aip 2 ) 
J y = n^latpi + ip\ aip 2 ) 

J d = ir: (^,0^2+ V4°Vi ) = ■ (9) 
They form the following subalgebra of the su(4)i current algebra: 



[4(x 


,J 3 M\ 


= 2iTr(e ijk J%(x) S(x 
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Sij S'(x 
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,J 3 M\ 
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(10) 



In terms of these currents the full Hamiltonian H — Hq + V H z + Hi can be expressed as follows: 

H = — dx J d {x) ■ J d (x) : + : J x (x) ■ J x {x) : + : J y {x) ■ J y {x) : + : J z (x) ■ J z (x) 

1 [°° 
H z = — / dx J z (x) 

^ Loo V ^ 

Hi = X d Jd(0) ■ S + XJ y (0) ■ S . 
The operator H z is the total (normalized) isospin z-component. [l|| 



(11) 



The electric current. We now turn to discuss in more detail the current across the quantum dot (we set the 
electric charge e = 1), and express it also in terms of the symmetry currents. The electric current J is given 
by the operator measuring the difference between the fermion density on, say, the second channel just before 
hitting the impurity and the fermion density on the same channel just after hitting it: 



J = lim (iplip 2 (x = -e) ~ i>\^i{x = 4) = ir~ Hm (J z (-e) - J z (e)) . 



(12) 



(we use here the "unfolded set-up" ) . Equivalently one can express the electric current as the rate of decrease of 
the charge on lead-2 (or increase on lead-1), 



J = -±N 2 = ± Nl 
dt dt 

= XJ x (0)-S. 



-i[H,H g ] 



(13) 
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It is easy to see that the two definitions coincide. We rewrite the first definition using "impurity conditions" - 
operator relation inherited from boundary conditions. Boundary conditions, in general, are part of the equations 
of motion and lead to operator relations valid on the full Hilbert space of a boundary quantum field theory. They 
are often derived from the action of the model in the same way as one derives the equations of motion. From 
our view point, after "unfolding" the Kondo model, we have a model with an impurity instead of a boundary. 
As with boundaries, impurities give rise to "impurity conditions" which are part of the equations of motion and 
are operator relations valid on the full Hilbert space. In operator language (which is more convenient for our 
purposes), the impurity condition associated to a local operator 0{x) can be written 



/ f — g poo \ poo 

lim I J dx + J dxj [H,0(x)} = J [H,(D{x)} 



Consider the impurity condition associated to the operator J z (x) with the model with Hamiltonian H ([2]) , 



l_im(J dx + J dxj [H,J z (x)}= J [H,J z (x)}. (14) 

On the left-hand side, only the free part Hq of the Hamiltonian is involved, because the operator J z {x) is never 
at the site x — (and H z commutes with J z (x)). Using the fact that with the free Hamiltonian Hq, J z (x) is a 
right-moving operator [Hq, Jz{x)\ = i-^J z (x), and using the asymptotic conditions J z (oo) = J z (—oo), we have 



lim / dx+ / dx) [H,J z (x)] = lim i(JJ-e) - JJe)) . 

^0+ VJ-oo Je J ^0+ 

On the other hand, on the right-hand side of HH), since the integration is on the full interval, the free part of 
the Hamiltonian does not contribute. Only the impurity term, at x = 0, contributes, and it gives 



/oo 
[H, J z (x)} = 2iir\J x (0) ■ S 
-oo 



as expected. 



Having now discussed the system and the various operators describing it we turn to discuss in more detail the 
nature of non-equilibrium in the system. 



III. EQUILIBRIUM VS. NON-EQUILIBRIUM 



Our model, a quantum impurity coupled to leads at different chemical potentials, describes a non-equilibrium 
situation - a current is flowing from one lead to another. What is the the precise meaning of this statement? 

In this section we show in what sense an out-of-equilibrium model differs from an equilibrium model. We begin 
by showing how the Keldysh formulation leads, when the system is in equilibrium, to the usual equilibrium 
density matrix description. We prove, in other words, the following: 



5 , (0,-oo)e- /3g °5 , (-oo,0) 
Tr [e-^o] 



-P H °Vexp(iJ° ip dtH\ a) (t) 



-I3H\ X 



Tr 



e-f>a°V exp (i dt fff (t)) ] Tr [e^*™] 



(15) 



as an equation to hold when evaluated inside traces with insertion of any number of local operators at fixed 
positions. 

A local operator is, by definition, an operator depending on the position x (in the sense that its commutator 
with the momentum operator is a derivative with respect to x) , such that its commutator with the hamiltonian 
density at position y is zero for x ^ y. Note that local charges, for instance conserved charges of the Hamiltonian, 
are integrals of local operators, and are not local operators themselves. Hence, the limit ([TB]) does not hold 
with insertion of local charges. This makes physical sense, since conserved charges are not expected to relax to 
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their equilibrium values. Technically, one must remember that the density matrix is an operator with infinitely 
many matrix elements, hence any limit applied to it cannot be expected to converge to an object having the 
same properties (or to converge at all) independently from which subset of matrix elements we are looking at. 

This derivation is important for what follows, so we present it in some detail. Recall that S(ti,t2), Eq. ([7]), is 
the evolution operator at zero voltage. We now establish some useful identities. In the interaction picture with 
respect to Hq we have, 

-iHot Oo -0Ho _ d „ VT , (a ^ TT<.°)(*\\ „-PH _ (■ [*° *" M ry-(O) , 




S(0,t )e~ lHoto e-P" = Vexp^i J dtH} u >(D)< -"••=, ' """P exp j / / ,UH,(I)\ (1(>) 

where in the interaction picture 

Hf\t) = e^Hje- 1 ^ = X d J d (-t) ■ S + XJ y (-t) ■ S . 

In (|15[) and in the last two expressions of (|16p . the symbol V indicates path-ordering in time: the operators are 
positioned from left to right with their time argument going from the lower integral limit to the upper integral 
limit. In the first occurrence in (fTI))) , integrals are ordered from on the left to to on the right. In the second, 
the integration contour is from — i(3 on the left to to — i(3 on the right. On the other hand, we have, 

JHotojSU n\ - -Do™ ( i [ A+ TfW 



S(to,0)<=Vexp\iJ dtHY\t)j . (17) 

The Keldysh evolution is then (multiplying (fTT|) with p6[) and dividing by the trace of this product), 

5(0,to)e-^ g °S(t ,0) 
Tr(e-^o) 





V exp 






Tr 


' e -l3H op exp fi jto-i/3 dt R {0) (t) 




)] 



\to\>/3 





Pexp (iSldtH?\t)) 




Tr 


'e-PBoVexp (i^ ip dtHf\t) 


) Ve^{i!ldtHf\t)) 



(18) 



The last equality is valid perturbatively if the integration from t — i/3 to to is negligible at every order in 
perturbation theory. 

Note that the last equality involves taking \to\ much greater than j3. At zero temperature, when (3 — > oo, this 
condition cannot hold, and since the correlation functions then may have algebraic decay with power —1 at 
large distances, our proofs below (at equilibrium and in the steady state) do not apply. Nevertheless, as will be 
seen, our two-loop perturbative results for the non-equilibrium current have finite zero-temperature limit; this 
will be discussed further in the last section. 

To show the last equality in (|18[) we evaluate the expectation value of a local operator (or product of any local 
operators at fixed positions) O, inserted at the right-hand side of the first equation of (|T5jl . Denoting by 

Tr (e-P Ho ■ ■ ■) 

Tr(eW ( 19 ) 
the averaging in the free theory at temperature /3 , we consider, 

[Ve^{iSX iP dtHf\t)) V^(iSldtHf\t)) O)) o 

(20) 



V exp (i dt H?^ (i)) V exp (i f° dt (t)) 

All correlation functions involved are correlation functions where the -ff/(i)'s are connected to O. Connected 
correlation functions are defined, in the usual way, by subtracting from correlation function appropriate products 
of expectation values. In Appendix [Al we recall their precise definition and main properties. Only connected 
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correlation functions occur, because in (|20p we divide by the correlation function of the operator where all the 
Hi(t)'s are involved. 

In Appendix [B] we show that correlation functions of the type 

((Hf (t + h)Hf ] (t + 1 2 ) ■ ■ ■ Hf ] (t + t n )O)) 

factorize, as t — > ±00, into 

((H^(t 1 )HP(t 2 )...HP(t n ))) Q ({O)) , 
with sub-leading asymptotic contributions vanishing exponentially for finite (3. 

The last step of (TT5)) follows from the factorization property. Due to this property, connected correlation 
functions of the type 

vanish exponentially whenever any subset of consecutive time variables {f,-, tj+i, . . . , tj} (corresponding to a 
subset of time-ordered operators Hj(t)'s) goes to negative infinity simultaneously. This implies that order by 
order in perturbation theory of l|20p . all integrands vanish exponentially in any large-time region, in particular 
in the segment to — which then factorizes and cancels between numerator and denominator. Hence, in the 
limit where t — > —00 the last step of (fT8]) is exact order by order in perturbation theory, and we have (fl~5j) . as 
claimed. See Figure 2. 

We wish to note that our argument in Appendix [B] relied on the fact that Hi{t) is a local, right-moving operator, 
and that it couples to the external degree of freedom S (the impurity) in an SU (2)-invariant way. For electronic 
degrees of freedom large time means large distance, and at large distances, correlation functions of local fields 
factorize. Combined with SU(2) invariance, this implies the factorization of correlation functions of Hj(t)'s at 
large |t|. 

The implications of the well-defined limit to//3 —> —00, and in particular of the factorization at large time 
separation of the correlation functions involved in the perturbative coefficients, are important. It was not 
necessary to invoke any external relaxation mechanism: the factorization signals a decoherence in time and 
suggests that Hq represents a good thermal bath, and this bath by itself provides such a mechanism. As 
in Caldeira-Leggett models, Hq can be seen as an infinity of free oscillators with an appropriate frequency 
distribution in order to represent a thermal bath. The loss of time-reversal symmetry associated with this 
relaxation mechanism occurs when taking the limit |to|//3 — > 00. 

The same derivation can be carried out for more general unitary conformal field theories perturbed by an 
interaction Hj at one point, or defined on a finite region of space. Inferring from our derivation, the interaction 
can be due to an external degree of freedom coupled to any linear combination of fields that factorize into 
right- and left-movers, and the coupling has to be invariant with respect to a symmetry group acting on the 
full configuration space. 

The steady state current. The derivation fails when out of equilibrium, V 7^ 0. The step that becomes 
incorrect, if we start with expression ([3]) for the steady-state average, is the shifting of the integration limits 

to — i/3 — > to. Indeed, the correlation functions involving VH z °\t) = VH Z are not suppressed at large negative 
times since H z is a conserved charge of the Hamiltonian Hq. Moreover, due to quantum fluctuations of the 
charge H z , made possible by the interaction Hj (that is, [Hi, H z ] ^ 0), connected correlation functions involving 
H z are not zero. These two conditions are at the origin of the appearance of a non-equilibrium situation. In 
physical terms, the first condition is that the bath represented by Hq does not provide a relaxation mechanism 
for reaching Boltzmann's distribution of states associated to the energy Ho + VH Z + Hi; the second condition is 
that nevertheless, H z is subject to quantum fluctuations and evolves with time. In this case then the Keldysh 
formulation does not reduce to an equilibrium description. 

Our analysis, however, has not yet established that a steady state occurs. We shall present below a full proof to 
this effect. To motivate the proof we begin with a physical argument, based on the first formulation described 
around Eqs. (HJ) and ([5]), by considering the respective ground states of Hq and H (instead of the associated 
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(A) 



to-ip 



-iP 




(C) 



to-iP 







-iP 



FIG. 2: (Color online) Deriving the equilibrium formulation starting from real time (Keldysh) formulation. (A) The 
Keldysh integration contour on the complex time plane; the directions of the arrows indicate the ordering of operators 
from left to right. (B) Adding the segment (to, to — i/3) is allowed by the factorization property (this step is not allowed 
when the system is out of equilibrium). (C) The equilibrium contour. 



thermal density matrices) and showing that they are "far" enough and that the evolution of H z is "slow" enough 
in the limit L — > oo so that a steady state is established. Under other circumstances, we might expect some 
oscillating behavior. 

That the ground state of H is far enough, and the evolution of H z is slow enough, can be made more precise 
in the following way. To begin with, consider the ground state |0) of Hq and the ground state \V) of the 
Hamiltonian Ho + VH Z . Later we shall consider the effect of the couplings A and A^. The ground state |V") can 
be obtained in the following way. Consider the operator 

U v =e i «f-<*> dxxJ * (x) (21) 
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(in order for it to be well defined, we assume an appropriate ultraviolet regularization of the operator J z {x)). 
This is a unitary operator, and its effect on Hq is: 

U_ v e~^ H °Uv — e~*^~ f dxxJ " ( . x ^ e ~ l3Ho e i ^ S dxxJ *( x ) 

— g— |3Hog- i% J dx xj z(x+i(3) J dxxj z (x) 

— g-/3-ff0g-«v / dx (x — i0)J z (x) f dxxj z (x) 

_ e -l3Ho e -f3VH z f dx 1 dx 2 (x 1 -i/3)x 2 [Jz(xi),J z {xz)] ^2) 

The last exponential factor is a real number, scaling with the system size L (it could be absorbed into the 
definition of U-v, a t the price of losing its unitarity). Hence, the operator J7_y takes Hq to Hq + VH Z , up to 
an additive number, and the ground state of Hq + VH Z can be obtained by 

\V) = U- V \0). (23) 



Computing the expectation values of H z (tracing over the two-dimensional impurity space) in these ground 
states we have, as L — ► oo, 

(0|fl,|0) = , (V\H Z \V) ~ -VL . (24) 

In particular, U—v has the effect, in the infinite- L limit, of changing the asymptotic conditions to J z (±oo) = 
— irV/2. We discussed the ground state of Hq + VH Z . However, the actual ground state of H yields corrections 
to the expectation values that are of higher order in the couplings with a finite limit as L — > oo, and our 
conclusion therefore apply to the full Hamiltonian. 

While the expectation values are infinitely distant in an infinite system, the rate of change of (H z ) as the 
interaction is switched on is finite since the operator Hi giving rise to the current is local. This will be seen 
explicitly in the perturbative calculations of the current J below (recall that J = —i[H, H z ]). 

Hence, as L — > oo, it would take more and more time to get from |0) to \V). Here we assume that the average of 
H z would decrease monotonically. This is expected for L large enough and elapsed time large enough, though 
not infinite. More precisely in the region ©, we expect the expectation varlue of H z to decrease steadily; this 
is the steady state. In other words, we expect a steady state to occur because H z scales with the length of the 
system, whereas its variation does not. For a finite L, it does not decrease monotonically at all times, and we 
might eventually see an oscillating behavior of period characterized by L. 

We proceed now to the main result of this section: we show that to all orders in perturbation theory the limit 
of very large negative times, to//3 — * — oo in ([5]), is well defined for any local operator, O, supported at a 
point or on a finite interval. This shows that there is indeed a steady state: the current operator J acquires a 
well-defined expectation value. 

Using the interaction picture with respect to Hq + VH Z , we can write the steady-state average of any operator 
O as 



(O)ss = lim — — — ^rrrTr 



V exp ^ Jdt H { p (tfj e-P Ho V exp [i Jdt H^ v) (t)j O 



o 

(25) 



with 

H [ p{t) = S H °+ VH ^ Hi e -<(ffo+vir.)t = e iVJr.t jff (o) (t)e -iVff,t _ (26) 
The operator Hf >(t) can be expressed in terms of "deformed" current-algebra operators. Consider 

jf\x) = e- iVH °*J d (x)e iVH * x 
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It is a simple matter to use the commutation relations (|10p in order to obtain 



Jd{x) 

cos(Vx)J x (x) + sm(Vx) J y (x) 
— sm(Vx)J x (x) + cos(Vx)J y (x) 



Using these operators, we have 



H 



{V) (t) = \djf\-t) ■ S + XJl v \-t) ■ S . 



(27) 
(28) 



The proof that the limit to//3 — ► — oo exists in (J5J) proceeds from arguments similar to those in the previous 
subsection. Let the operator O be supported on a finite interval. Moving the operator e~^ H ° to the left inside 
the trace in (|2"5T). we have 



where 



H<P(t) 



lim 

t //3^-oo 



Tr 


e-0 ff ° V exp (i f°~ iP dtH[ V) (*)) 


Vexp{iJ t o dtH^(t)) 


O 


Tr 


e -/3ffc V exp (i dt H { p (t) 




)\ 



(29) 



?-0VH z tt(V) 



f3VHz 



H\ '{t)e 

\ d J d (-t) -S + Asin(V \t + ip))J x (-t) ■ S + \cos(V(t + ij3))J y (-t) ■ S 



The exact form of Hj is actually not important; note only that it is a linear combination of local operators evolved 

in interaction-picture time. Again using the fact that only connected correlation functions (where Hj (t) and 

Hj (t) are connected to O) occur order by order in perturbation theory, and the fact that correlation functions 

involving Hj V \t) and Hj(t) factorize at large times t, one can see that all integrals are convergent in the 
limit to/ p — > — oo order by order in perturbation theory. 

Physically, this means that the bath represented by Hq provides the same mechanism for the steady state to 
occur as the mechanism it provides for the system to reach equilibrium in the case V = 0. 



We now cast the expression for the steady state averages in another suggestive form and derive the alternative 
formulation, expressed in ([5]), with the steady state obtained by coupling the dot (i.e. turning on the couplings 
A, Xd) to leads initially equilibrated at temperature T and at potential difference —V. 

Observe that the operators with superscript (V) form the same current algebra, Eq. (|10p . as those without 
superscript since they are obtained by the unitary transformation Uy (|2T1) : 



jy\ X ) 



U- V Jd(x)Uv 
U- V Jx(x)U v 

U-vJy{x)Uv 



4 v \x) = U-vJz(x)U v = J z {x) + ^ 
Hence, the steady-state average of an operator O can be written 



(0) s 



Tr 



lim 



Vexp(if*°dtH<?\t)) Uve-P Ho U- V rexp(if? o dtH<?\t)) U v OU- 



Tr [e-f JH "} 



(30) 



(31) 



where we recall that H^p (t) is the operator evolved with Hq only. Recalling the transformation of Hq under 
U v (122), we find 
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where 

{ - v) = UyOU-v ■ (33) 

In (|3"!2)) . it was necessary to include the factor e^ VH " inside the trace in the denominator of the right-hand side. 
For a system on a finite interval, the inclusion of this factor has the effect of cancelling the constant term that 
appears in (|22p . Then, the limit of infinite interval is well defined. 

Since for the current operator (fT2"|) we have J^~ v > = J , this shows the equivalence, for the steady state current, 
between the formulation §5§ and the formulation ([8]). In general, we will denote the steady state average in the 
latter formulation by 

{0)ss ' = W^-co Tr [r«W] ' (34) 

That is, 

(0) ss , = (U-vOU v ) ss ■ (35) 



Below, we carry out some formal manipulations which are justified only if we can establish a more stringent 
convergence property as that used above. We need to establish that the following expression: 



(0} ss i = lim lim 

t /@—*—oo t' o /0—t—oo 



Tr [5(0, f ) e -K H °- VH J S(t ,0) O^] 



Tr[5(0,t' o ) 



-13(H -VH Z ) 



S(to,0)] 



(36) 



with the limits on to and on t' taken independently, will yield a result independent of the order the limits were 
taken. Note that we have included factors S(t o ,0) and S(0,t' Q ) in the denominator. They assure convergence 
and cancel by cyclicity of the trace if the limit exists. To prove convergence in (|36l) . we consider 



(0) ss — lim lim 



Tr 


Pexp Uj*°dtH( V) (tj) 




O 


Tr 


"Pexp 




) e-^Pexp^/^ dtH { p(t) 


)] 



(37) 



Indeed, the same arguments we used to establish the connectedness and factorization allow us to take, for 
instance, first the limit with |f | large, then the limit with \to\ large, or vice versa. The result is unique. Using 
the operator Uy in a manner similar to the one above, it is a simple matter to obtain (|36[) from (|37[) and the 
result then is the same as the one obtained from the formulation (l34l). 



Alternative description of the Steady State. What replaces the density matrix e~° H description of 
equilibrium? We could translate the proof establishing equilibrium when (V = 0) to the proof establishing 
steady state when (V ^ 0) by means of the current algebra of symmetries. By similar means we shall show 
that a new operator will play for the system in its steady state a similar role to the one played by the density 
matrix in equilibrium. Such a steady-state density matrix can be obtained from simple manipulations, now that 
we have established the convergence of the integrals. 

Consider the formulation (|3l))) of the steady-state problem, with O in (f3"3"|) an operator supported on a finite 
interval in the theory Ho. Bringing e~P H " completely to the left, the right-hand side of (f3"6")l can be written as 
follows: 



lim lim 

to/0— *— oo t'o/P-^-oo 



Tr 




e^Pexp (ifcdtH?\t)\ 




Tr 




(*)) 


efWB'Vexp (iJ^dtHf* 


it)) 





Since we showed that the limits can be taken independently, we can shift t' — i/3 to t' Q both in the numerator 
and in the denominator, without shifting io, with vanishing error in the limit. We can then take t — to and 
keep only one limit symbol. Inserting 



1 ='P<-K\ i [i 1^ dtHPit^j "Pexp (ij^ ° dtHf\t) 
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just before the operator e^ VHz we have, 



Tr 

lim 



Defining the operator: 



-0H\v =o V exp U jj° dt ijj 0) (t)) e? VH * V exp (i /° ffj 0) (t) 



F= Urn £(0,i o )ffz£(i o ,0) (38) 



allows us to write the steady-state average of a local operator as 

_ Tr [ e -^lv=o e^O] 
{ ° )ss ' = Tr [e-^»lv=o e /3vr] ■ ( 39 ) 

Note that the limit (f3"5| cannot be expected to exist as an operator (recall that we are dealing with operators 
with infinitely many matrix elements), but only when inserted into appropriate traces (or only when appropriate 
matrix elements are considered) . More precisely, we have only proven that ([38)) is a well-defined operator when 
it is evaluated in expressions like (l39|) . and that the result is the steady-state average of the local operators 
inserted. This is a statement solely about a small part of the matrix elements of the operator ([38]) . The meaning 
of Eq. (|3"9")l is that one must first evaluate the traces and their ratio with the expression {55J at finite to/ ft, then 
take the limit indicated in ([38]) on the result. The properties of Y as an operator acting in a Hilbert space will 
be discussed elsewhere. 

Observe, however, that in all situations where the operator Y is well defined (that is, when we consider the 
appropriate matrix elements), then it is a conserved charge. Indeed, when it is well defined, then the limit t /f3 — > 
— oo of 5(0, t )H z S(t , 0) (or of any function of this operator) must exist. Since 5(ti, t 2 ) = 5(ti + dt, t 2 + dt), 
we have 4z lim t0 / (3 ^_ oo S(t, to)H z S(t 0: t) = 0, hence [H |y=0j Y] — 0- Then we can finally write 

Tr \ e -t3{H\ v = -VY)Q\ 

= TV[e-*g|v-o-v y)] J • ( 4 °) 
That is, averages in the steady-state can be obtained by tracing with an appropriate density matrix [20j. 

What difference is there between the equilibrium and the steady state? Consider a quantum mechanical system 
described by a Hamiltonian H. Put the system at equilibrium with a bath where there can be exchange of heat 
and of any quantity Q that is conserved by the dynamics H. The average of observables is then described by 
the density matrix e~ l3 ^ H+ ' 1 ^ where [i is the chemical potential associated to Q: the energy brought to the 
system by increasing Q by one unit. In expression (|4TJ|) . the average of a local operator O in the steady state is a 
trace with a density matrix of exactly the same form. The main difference is that the operator Y is a non-local 
conserved charge. A local conserved charge can be written as an integral over space of a local operator of the 
theory H\v=o plus a local operator at the impurity site, with possible non-trivial impurity-space components. 
The operator Y (I3"51) cannot be written in that way. To see this, we can write it as follows: 

Y = H Z + f dtj(t) (41) 



where 

J(t) = S(0,t)JS(t,0) (42) 

is the time-evolved current J with respect to the theory H\v=o- Then it is simple to write it as an integral of 
a charge density: 

/oo 
dxj tot (x,0) (43) 
-OO 

with 

j tot (x, t) = -J z (x, t) + S(x) [ dt' J(t') . (44) 

7T ./_oo 
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The charge density has a local bulk part, but the term at the impurity is not a local field of the theory H\v=o- 
it is the time integral of the current, and the current is not the time derivative of a local field. 

The non-locality of Y is the main difference between the description of a steady state and of an equilibrium 
state. In the formulation (|4"D")) . only a restricted set of operators O have well-defined average: those that have 
stationary expectation values. All local operators are of this type, but, for instance, it is simple to see that the 
operator H z does not have a well-defined steady-state value. 

Note also that the operator Y gives in principle a description of the asymptotic state that one can use in order 
to describe the steady state: quantities in the steady state can be evaluated as averages in an appropriate 
asymptotic state. Further analysis in this direction will be presented in our future works. 

We wish to remark that sometime ago Hershfield @ has considered steady-state flow and has argued that under 
some assumptions concerning the relaxation of correlation function an expression (|40[) would govern the steady- 
state current. He gave then implicit equations to determine Y. It appears to us that our explicit expressions 
for the operator Y satisfies his implicit equations, and should hence correspond to the same operator (although 
we have not thoroughly ascertained the confluence of the two approaches). We must stress, however, that no 
assumptions were made in our derivation. 



IV. RG-IMPROVED REAL-TIME PERTURBATION THEORY 



The perturbative expansion. We now turn to real-time perturbation theory for the current 1(5)). We take the 
formulation where the system is initially brought to equilibrium with a nonzero bias voltage, then disconnected 
from the external bath before the voltage is turned off and the interaction is turned on. It will be convenient 
to consider adiabatically turning on the interaction in the infinite past: we introduce a large-time exponential 
cutoff, e vt Hj, with 77 a positive scale with dimension of energy, and take the limit to/ (3 — ► — 00 in ([5]). The 
quantity rj/3 will be sent to at the end of the calculations. This means, physically, that the two leads are 
slowly brought towards the dot after the voltage has been turned off. Our proof that there are no divergencies 
as to/P — > —00 in the previous section shows that there are no divergencies as r\(3 — > 0. The current can then 
be written 

{J)ss = Jr 0+ Tr [e-Mo-VH.)] ( 45 ) 



with mj 



More precisely, 



S n (h,t 2 ) = V exp [' i(H + e*Hi) dt . (46) 
Jt-i 



00 M M rO 

(J) ss = lim Vi fc dhe^ dt 2 e^ ■ ■ ■ / dt fc e** fc (([ff/(*i), [Hj(t 2 ), • • • , [#/(**), .7] • ■ ]}))v (47) 



k=0 

where 



Tr r e -P{Ho-VH z ) . . .] 
((■■■))V= r^rp^vH^ ■ ( 4 8) 



The integrals in this expansion are plagued with ultraviolet divergencies which we have to regularize. This can be 
done in several ways. For our purposes, it will be convenient to modify the operators Jd(x), J x (x), Jy(x) in order 
to render their correlation functions regular at coinciding points. More precisely, we choose the regularization 
scheme where all operators (in the Hamiltonian and in correlation functions) at the impurity site are regularized, 
whereas all operators away from it are unaffected. Since the interaction is only at the impurity site, this is 



17 



enough to regularize the theory. Define the momentum space (mode) operators Jd(p), Jx{p), Jyip) and J zip) 
in the following way: 



J d (x) = dp.J d (p)e lpx , J x (x) = dpJ x ( P y px 

J — oo J -co 

pOO p CO 

J y (x)= d P j v ( P y px , J z (x)= d P j z { P y px 



The mode operators satisfy the following set of commutation relations: 

[Jdiv), J j M\ = 1 e V* J d(P + <l)+P % S(p + q) 



[J*(j>), 4(q)} = i e ijk J k d (j> + q) +p5 l:j 5(p + q) 

[J l y{p), Jib)] = * eijk Jd(p + q)+p % S(p + q) 

[Jz(p),Jz{q)} = pS(p + q) 

[Jd(p),4b)] = i^kJ k x {P + l) (49) 

[Jd(p)>4(a)] = Jy(p + q) 

[Jl(p),JM\ = o 

[JliplJiiq)] = Jz(p + q) 

[Jy(p),Jz{q)} = iJ x (p + q) 

[Jz(p),J x (q)} = iJ^ + q). 
We then introduce the regularized operators 

/CO pOO 
d P R A (p)j d ( P y px , (f x ) A (x) = / d P R A (p)j x ( P y px , 
CO J —oo 

CO 



(j v )a(x) = / dpR h {p)j y (py px 

where R A {p) is a function that vanishes as \p\ — > oo, and whose complex conjugate satisfies R A (p)* — R A (—p) 
in order to preserve hermiticity of the regularized operators. The function R A (p) can be chosen in many ways, 
and the choice is a matter of convenience. We will choose a gaussian regularization, 

R A (p) = e"P 2 /(2A 2 ) _ (50) 

The parameter A plays the role of an effective band width (we do not denote it D since it is certainly not exactly 
the band width). The universal part of the limit A 3> V,T is the same as that of the limit of large band width. 

The time integrals in (|47p can now be traded to momentum integrals: 

dpi_R A {pi) f dp 2 R A (p 2 ) f dp k R A (p k ) 



, n J Pi + ill J Pi + P2 + 2iv J 



Pi + P2 + ■ ■ ■ + Pk + ki-q 
f dsR A (s)(([H I (p 1 ),[H I ( P2 ) 1 ---,[H I ( Pk ),J(s)}---}})) v (51) 
where 

Hi(p) = X d J d (p)-S + 2J y (p)-S (52) 

and 

J{p) = J X {P)-S. (53) 

The parameter rj can be set to zero with the requirement that the momentum integrals be taken on a line 
parallel to the real axis in the p-plane with a slight positive imaginary part: 

dpiR A (pi) f dp 2 R A (p 2 ) f dp k R A {p k ) 



,jx = V^(_X) fc [ dl>1 Ra (P^ f d P2 R A(P2) f 

tr'o J+ pi J+ pi +P2 h 



Pl+P2-\ h p k 

dqR A (q) (([H I (p 1 ),[H I (p 2 ),---,[H I (p k )J(q)]---]})) v . (54) 
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Note that this can be done since Ra(p) does not have a singularity at p — or anywhere on the real p-line. No 
singularity at pi + ■ ■ ■ + pj =0 (for j = 1, . . . , k) can occur in the averages ((■ ■ -}}v hi the expression above, 
since we know that there are no divergencies as f][3 — > 0. 

In fact, for explicit calculations, it will be more convenient to integrate over the real line in momentum space, 
and to use the formalism of principal value integrals. For this purpose, recall that 

dp . f{p \ = [dp f ( ~ 1)n ^("-i)(p + g )+p 1 - \ f{p) (55) 
+ (p + q) n J l(n-l)! ^ -(p + q) n J 

where P says that we must take the principal value of the integral: 

J dp ^P-^ — ~~^y^j f(p) — finite part of (^J +^ ^ dp ^ ^ in power expansion as e — > . (56) 

Explicit one-loop calculations. The trace over the impurity space of the operators 

[H I ( Pl ),[H I (p 2 ),---,[Hi(p k ),J(q)}---}} 

can be obtained by using the following general formula, valid for any vector operators A and B that commute 
with S: 



[A-S,B- §\ = t {A\B*}e ijh S k + \[A\B l ] ( .-,7i 



where {•, •} is the anti-commutator. Using the commutation relations (|49[) for evaluating the commutators, one 
is left then only with multiple anti-commutators of the mode operators. To one loop we will need the following 
commutators: 

Ci = [J y {q)-S,J x (k)-S] 

C 2 = [J d (p)-S,[J y (q)-S,J x (k)-S]} 

C 3 = [J y (p) ■ S, [J d (q) ■ S, J x (k) ■ S}} . (58) 
Under trace over the impurity space, they give 

Tri mpur ity(Ci) = -3iJ z (q + k) 

Tr impurlty (C 2 ) = 2{Ji(q),Ji(p + k)}-2{Ji{ P + q),4(k)} 

Tr impurity (C 3 ) = -2{.r y {p + q),Jl{k)} . (59) 

The traces ((• • -))v over the bulk-CFT Hilbert space of these anti-commutators can be calculated without the 
need to construct the Hilbert space for Hq — VH Z , but only by using the following exchange relations: 

J d {p)e- p{Ho - VH ' ] = e-^- VH ^ J d {p)e~^ 
J z (p)e-^ H °- VH ^ = e-Wo-VH,) j z{p)e -0 P 

J + (p)e-^ H °- VH ^ = e-^ H °- VH >Kj + (p)e~^+ v ^ 

J_(p) e -0(Ho-VH z ) = e -p {Ho -VH*) J_( p)e -P( P -V) (6Q) 

where we use the following linear combinations: 

J + (jp) = \(Up)+iJy{p)) > j-{p) = \(j x {p)-iJy{p)) • (61) 
The exchange relations imply for any operator O 

((OJ d (p))) v = ((J d {p)0))v 

({OJ z (p))) v = e-e* ((J x (p)0))v 

((OJ + (p))) v = e-^+v) ((J + (p)0))v 

{{O.L{p)))v = e-^- v H(,L(p)0)) v . (62) 
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Using the commutator formulas (l4"9")) , using these relations and using the traces of single operators, all traces 
can be calculated. Traces of single operators can be calculated from the fact that ((Jd,x, y ))o = ((Jz))o = 0, and 
from {{0))v = ((U-vOUv))o- In terms of modes, the operator U-y has the representation 



U- V = e 2VJ '^ 

where we formally define J' z (p) = dJ z (p) / dp. The traces of single operators are then given by 

{{Jd,w, v (p)))v = 0, ((J z (p))) v = V5(p) . 
The trace of the only type of anti-commutator appearing in (|59[) is then given by 



«K(p), 4(<l)}))v = -r;F(p)6(p + q)6 i 



where 



From ([54")) and using (|55|) . the integrals to be calculated at one loop are 

r 



(63) 



(64) 



(65) 



(66) 



h 
I2.3 



dqR A (q) -in6(q)+P- 



dkR A (k) «Ci)> v 



J d P R A (p) (-iirSfr) + J J dqR A {q) \-iir6(p + q) + P^-j- 
dkR A (k) ((C 2 , 3 ))v ■ 



The current will then be given by 



(J)ss = A 2 {h + \ d (h + I 3 ) + 0(X% A 2 )) 



(67) 



(68) 



Since all integrals are real, it is clear that in the expressions (f6"7| , only terms with an odd number of delta 
functions in the momentum variables give non-zero contributions. It is a simple matter then to obtain 



Ii = 
h=h = 



3ttV 
3irV 



dpL-f(p) Ra(V P ) 2 
P 



where p is now a dimensionless momentum variable, and 



/(p) = (p + l) 



and 



1 + er w{ -P +1 "> 
1 _ e -tu(p+l) 



w =0V 



(P-1) 



1 + e-""^- 1 ) 

1 — e -w(p-l) 



(69) 
(70) 

(71) 
(72) 



Using the symmetry f(—p) — —f(j>), the integral Ii can be calculated in the following way, keeping only the 
divergent and finite parts as A — > oo: 



3ttU 
~Y~ 
3ttU 
~Y~ 

3ttV 



Jo P 



dp 



o P 
P(w) + In 



(/(p) - 2 + 2e-* ) - 2 



* ( e -p 2 _ i) e - P V 2 /A 2 



o P 



(73) 
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where the symbol ~ means that equality is valid only for the finite and divergent parts as A — > oo and where 

P + 1 p—1 _ 



This is a well-defined function for all w with positive real part, and is easy to evaluate numerically. 



(74) 



It will be useful to have the asymptotic behavior of the function P(w), at large and small w. This is evaluated 
in Appendix [D] The asymptotic expansion at large w is given by 



P(w) - 1 + 



7 



-w 



as w — > oo , 



(75) 



whereas the expansion at small w is 



P(w) = ln(iu) + ^- + 7 - ln(27r) + 0{w) as w -> 



(76) 



Summarizing, the zero- and one-loop divergent and finite contributions to the current are (with T = 1//3) 



p+1 



with P(io) given by the integrals in (74|), P(w) = J °° ^ (- e „ (p+1) _. . ,. 
combination P{V/T) + ln(A/V) appears at one loop. This combination has the limits 



(77) 

Note that only the 



1 + 1 



-|l,,|T)-^- 7 + ln(to) 



as T < V < A 
as V < T < A - 



(78) 



The one-loop calculation was also performed in [l7T | (although the analysis did not go in as much detail as ours), 
and it can be verified that their results agree with ours. 

Results for the two-loop calculations. Two-loop integrals come from the commutators (here p, q, r and k 
are all momentum variables): 

d = [J d (p) ■ S,[J d (q) ■ S,{J y (r) ■ S,J x (k) ■ §]}} 
c 5 = [Up) ■ S,[J y (q) ■ S,[J d (r) ■ S,J x (k) ■ 5]]] 
C 6 = [J y ( P )-S,[J d (q)-S,[J d (r)-S,J x (k)-S]}} 

C 7 = [J y (p)-S,[J y (q)-S,[J y (r)-S,J x (k)-S]}} (79) 

appearing inside the traces in the integrands of (|54"]) . Calculating the traces and using (|55[) . we can find the 
corresponding set of two-loop integrals, l4,I^,Ie and I7 (written in Appendix [E| . These integrals enter the 
current as 



{J) ss - A 2 (h + X d (I 2 + h) + X 2 d (h + h + h) + + 0(X% X 2 X d )) 
Explicit calculation of these integrals lead to the following divergent parts as A — > 00: 



(80) 



h + h 



h 



3ttV 
3ttV 



IOP(io) +5 In 



- 1 



In 



finite 



2P( W )+ln(-) -1 



A 



ln[ - 

V 



finite 



(81) 



The finite contributions are much more complicated, and are reported in Appendix [FJ Let us only notice that 
the limit T/V — > of these contributions is finite, as in the one-loop results. This indicates that the voltage 
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V plays the role of a good infrared cutoff for the perturbative calculation of the current to two loops: the 
temperature may be set to zero without divergences to this order. 

Summarizing, the zero-, one- and two-loop finite and divergent contributions to the current are, 



(J),, = <j \ I 2A d (p(w) • In f - 



+A 



"10/ V') -h '.In ( ^ ] - I 



2P(>i + ln(~) - J 




y / 

+\ 2 \ d [h + h + h] finite + A [l r ] finite 

(82) 

where [h + h + ^finite and [^7] finite are given in Appendix [F] 

The Renormalization Group equation (The Callan-Symanzik equation). In systems at equilibrium, 
Wilson's renormalization group ideas allow us to understand how physical quantities can have universal forms 
(independent of the precise form of the interactions at the microscopic level) when all physical energy scales 
(temperature, voltage, etc.) are much lower than the microscopic scales (band width, inverse lattice spacing, 
etc.). Out of equilibrium, it is not obvious that Wilson's renormalization group ideas still apply. 

A different, but equivalent, way to look at universality is to study the limit where the cut-off A is increased and 
sent to infinity, while the couplings dependence on the cut-off is again governed by an RG equation, valid at 
very large cut-off. If such a limit exists then all quantities tend to their universal form. 

We will argue that the steady-state average of the current operator (or of any operator having a well-defined 
average in the steady state) satisfies the Callan-Symanzik equation with the same beta function and anomalous 
dimension as they occur in any average evaluated at equilibrium. More precisely, we will argue that 

+ /MA,A rf )— + /?A d (A,A d )— ){J) SS A =°°0 (83) 
\,\d dX dXd ) 

where /3a and 0\ d are the beta functions of the anisotropic two-channel Kondo model. Note that the anomalous 
dimension term does not occur: the current operator J has zero anomalous dimension (this is natural from a 
physical perspective, as the current is a physical object which should not change with a change of scales; we 
will verify this explicitly to one loop, and indirectly to two loops, below). 

We should note that the Callan-Symanzik equations with one-loop beta functions and zero anomalous dimension 
was written in (l7j for the steady-state average of the current, from physical arguments (22I]. Here we present a 
quantum field theoretic argument that applies to all orders (and all matrix elements) , and in the next section we 
explicitly verify this argument and calculate the beta functions and anomalous dimension in a universal fashion 
(so that it automatically applies to the steady state) to one-loop order. 

The Callan-Symanzik equation embodies Wilson's renormalization group ideas: it tells us how a change of cutoff 
A (for instance, the band width) can be compensated by a change of few relevant coupling constant, as long as 
all physical energy scales are much lower than A. Solving the Callan-Symanzik equations (this is done below 
for the current in the steady state) allows us to describe the low energy behavior of the steady-state current in 
terms of the ratios V/Tk and T/Tjc, where Tk is an integration constant, as well as of one extra parameter 
(invariant under the RG flow) characterizing the asymmetry between the couplings \ d and A; we will denote 
this parameter by C. The integration constant Tk and the extra parameter C characterize the quantum field 
theory; when they are fixed, all averages can be evaluated unambiguously. These parameters are not universal: 
different microscopic theories have low-energy behaviors described by different values for them. Up to these 
non-universal quantities, the quantum field theory description is universal, independent of the precise choice 
of the cutoff procedure (precise structure of the band, for instance) and of irrelevant couplings (interactions 
that give vanishing contributions at low energies). The integration constant Tk is the Kondo temperature: 
the temperature above which the "Kondo cloud" gets destroyed by the thermal energy. It is related to the 
microscopic values of the couplings A and Xd (the values when A is chosen to be of the order of the real band 
width), and it decreases if the couplings are decreased. At zero couplings, the Kondo temperature is zero and 
all scales of the low-energy physics disappear: this is a quantum critical point. The quantum field theory with 
finite ratios V/Tk and T/Tk describes the situation Tk <C A: the couplings are sent to zero at the same time 
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as the voltage and temperature are made much smaller than A. This is the scaling limit, describing the region 
around the quantum critical point A = A<j = 0. In particular, one finds that the perturbative expansion is valid 
in the region Tk < V or T K < T. 

We will see that in the steady state, it will be more convenient to introduce a scale Mk characterizing both the 
effects of the thermal energy and of the electric potential driving the current on the Kondo cloud. We will then 
be able to compare these effects, using a comparison parameter that is exact in one-loop perturbation theory. 

It is natural that the Callan-Symanzik equation is still valid in the steady state, since the steady state can 
be understood as an appropriate asymptotic state, characterized by a scale V, and since averages of operators 
in asymptotic states satisfy the Callan-Symanzik equation. In particular, V should flow trivially with the 
renormalization group. However, it is instructive to see explicitly how this works in real-time perturbation 
theory. 

We will argue that Eq. (|83[) holds simply from the fact that the Callan-Symanzik equation is satisfied for any 
average in equilibrium. The main observation is the following. Consider the Hilbert space H-v associated to 
the Hamiltonian Hq — VH Z . It is not formed of vectors that are in the Hilbert space Ho associated to the 
Hamiltonian Hq. In particular, its ground state | — V) can be formally defined as 

Uv\0) , 

where |0) is the ground state of Ho- This definition indeed makes sense for any finite length of the system, and 
in fact allows to calculate matrix elements of any operators also at infinite length, but it does not give a vector 
in Ho at infinite length. Nevertheless, the mode operators associated to current algebra operators still have a 
well-defined action on H-v ■ In order to see this, it is convenient to construct the Hilbert space Ho by defining 
a vacuum state |0) satisfying J +i _ i£ j(p)|0) = and J z (p)\0) = for all p > 0, and by constructing other states 
of the Hilbert space by acting with J + ._ j( z(j?) and J z {p) at p < 0. Then, it is a simple matter to see, using (163)) . 
that 

J+{p) \ - V) = if and only if p > -V 
J-{p)\ - V) = if and only if p > V 
■Jd(p)\ - V) = if and only if p > 

Jz(p)\ - V) = if and only if p > . (84) 

The main observation is that any normal ordering operation valid on Ho is still a good normal ordering operation 
on H-v ■ Indeed, a state of the form 

j^( Pl ) ■ ■ ■ r + «( Pa ) j^(pi) • • • J- b (p' b ) ^(Pi) • • • Jl e (p") UpT) ■ ■ ■ Mp7) I - v) 

gives zero whenever 

pi + . . . + p a +p[ + . . . + p' b + p" + . . . + p" + p'" + . . . + p'H + (a - b)V > . 
Since (a — b)V is a much smaller than A for any finite a and b, a normal ordering of the type 

■J(pi)J(P2)~-J(pk)'-=J(pi 1 )J(pi a )~-J(Pi k ) with p h <p i2 < ■ ■ ■ <p lk , (85) 

where the indices i m 's are all different and drawn in appropriate fashion from the set {1, 2, . . . , k} as to make 
the set of inequalities for the Pi m 's valid, is still a good normal ordering on H-v 

This observation is enough in order to see that the Callan-Symanzik equation is still valid in the steady state. 
Indeed, the Callan-Symanzik equation is really an equation for operators, rather than just for particular averages. 
Recalling the regularized real-time perturbation theory (|54l) . consider the operator 



Ja = £(-1)"'/ d Pl R A ( Pl ) (-iir6{px)+P^J J dp 2 R A (p 2 ) (-iTrSfa +p 2 ) 
• • • x J dp k Rxipk) ^-«7r<5(pi +p 2 H \-pk) + P 



P2 

1 



'Pi + P2 H VPk 



J dsR A (s) [H I (p 1 ),[H I ( P2 ),---,[H I (p k ),J(s)]---]] . (86) 
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It gives the interacting current operator with respect to the free Hilbert space Tto (at equilibrium) or H-v (hi 
the steady state) . The validity of the Callan-Symanzik equations at equilibrium means that all matrix elements 
of J a in the Hilbert space Ho satisfy the Callan-Symanzik equations. This can be written: 




A-^ 



+ /? A (A,A d ) — + p Xd (\,\ d )— ) = on H . (87) 



Indeed, such matrix elements correspond, in perturbation theory, to matrix elements of the current operator J 
in the basis of eigenstates of the full Hamiltonian H\v=o (all of which obey the same Callan-Symanzik equation). 
We give an argument for the validity of ([57]) in Appendix [UJ In ([57]) . the limit A — * oo should be performed 
after a matrix element has been calculated, holding the momenta associated to this matrix element fixed. We 
now extract the divergent and finite part of the operator J as A — + oo: 

Ja-E^' (A) (88) 

where Oj 's are operators built out of the mode operators for the current algebra along with possible impurity- 
space operators (the normal ordering is defined in (|85|) ). This is always possible to do by rewriting Sa in terms 
of normal-ordered operators and evaluating the coefficients as A — > oo. The Callan-Symanzik equation (|87[l 
then states, 

:O j+1 : + (p x (\,\ d )-^+p Xd (\,\d)^j :O j : = (j = 0, 1, . . . , oo) . (89) 

Since the normal ordering operation (|85[) is also valid on the Hilbert space Ti-v, it is clear that the expression 
([88]) also gives the divergent and finite part of the current operator on Ti-v, so that the recursion relation ([89)) 
among operators : Oj : implies the Callan-Symanzik equation (|83[) also holds for the steady-state average of the 
current. Of course, the same is true for any operator that has a well-defined average in the steady state. 

All this will be explicitly verified for the current operator to one loop in the following sub-section. 

Density-matrix-independent calculation of the beta functions and of the anomalous dimension of 
the current to one loop. Let us write Sa using time variables instead of momentum variables: 

JA = J2 ik dh dt 2 --- dt4(H^)A(t 1 )A(H^)A(t2),---d(Hf ,) )A(tk),(J (0) )Am---}} ■ (90) 

f,_Q J —oo Jt\ Jtk-x 

Here 

{Hf ] ) A {t) = (\{J y ) K {-t) + X d (J d ) A (-t)) ■ S (91) 



and 



(J i0) h(t) = X(J x ) A (-t) ■ S . (92) 



In the regularization scheme that we consider, characterized by the regulator Ra{p) (|50p . it is a simple matter 
to observe that 

a Iv ( ^ 0))aW = -^ H i 0) y^ > a Jv (i7(0))a(<) = - a^ (0) ™ (93) 

where primes mean time derivatives. Consider the first few terms of (|90[) : 



J a = {3^)a + / dt [i(H?%{t), ( J w )a(0)] + 

J — oo 

dt! [ dt 2 [*(HW) A (ti), [i(HP) A (t 2 ), (.7 (0) )a(0)]] + . . . (94) 



— OO 
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Using integration by parts, we find 



OA A 2 



*(ff| O) )A(0), (<7 (0) )a(0)] + / dt [i(#} 0) ) A (i), (J^)l(O)} + 

J — CO 

^[[ l ( J ff| 0) )^W, i (iji 0) ) A )(t)],( l 7(°)) A (0)] + 
dt[(Hf ) ) A (t),[(H ( I 0) y A (0),(jW) A (0)}} + 

) 

dh f dt 2 [i(ifj 0) )A(«i), [^fkfe), (.7 (0) )a(0)]] + ...V (95) 



— oo 



We want to evaluate all this at A — > oo in order to find the beta function using the equation (|57|) . Inside the 
parenthesis, we need only keep the contributions of order A 2 In 7 (A) for non- negative integers j. The very first 
operator of course does not contribute, but all others do. These leading contributions have to be compared with 
the leading contributions of 



(96) 



where we wrote only the terms contributing to the one-loop order. The beta functions appearing there can 
be obtained by requiring that when the derivatives with respect to the couplings are applied to the operator 
i(Hj°^) A (t), they give the operator [i(Hj)' A (t),i(Hj) A )(t)] (appearing on the third line of (pa)) ) in the limit 
A — > oo. The main contribution in this limit of this operator can be obtained from 

±{i(HP)' A (t),i(HP) A )(t)] = -i ((A3 + A 2 )(J d )^ A (-<) + 2\\ d (J y ) v - 2A (-t)) -S+j^-.O: (97) 

where O contains products of current algebra operators at the same point x = — t; the explicit form of O is not 
important here. Equating the leading behavior of this operator as A — > oo with that of 



gives 

/?A d =-(\ 2 d + \ 2 ) + 0(\l\\\ 2 d \ 2 ) , (3 x = -2XX d + 0(XXlX 3 ) . (98) 

Note that we had to take the limit A — > oo of the operator [i(Hj )' A (t), i(Hj )h)(t)] rather than of the integral 
where it is involved on the third line of (f95|) . since the beta function does not depend on the particular average 
that we are calculating. The sub-leading operators in (|97j) may give contributions to this integral, but these 
are two-loop contributions to the anomalous dimension of the operator J (of course, since this operator should 
have zero anomalous dimensions, all such contributions should cancel out). 

The total one-loop contributions to the anomalous dimension of J can be verified to be zero by checking that 
when the derivatives with respect to the couplings in (|96[) are applied to the operator {J^^)a(0) (the first term 
inside the parenthesis), they give the two terms appearing on the second line of (|95[) in the limit A — ► oo. In 
this limit, the second term on the second line of (|95p can be written 

" dt[i(H^) A (t),j"]r. f dt[i{Hf ) )l{t),J] = \i{H? ) )' A {a),J}. (99) 

-oo J —oo 

That is, it is equal to the first term. Together, their leading behavior at A — > oo can be obtained from 

^[i(Hf)' A (0), (J^MO)} = -2X^)^(0) + ^:0: (100) 
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where O is of the same form as O. But since 



we immediately conclude that the anomalous dimension of the current is zero to one loop. 



It is a simple matter to verify that the one-loop steady-state current (|68|) with the values (|69|) and (|T3[) indeed 
satisfies the Callan-Symanzik equation 



A,A rf 



(101) 



Note that our derivation of the Callan-Symanzik equation to one-loop did not use the particular initial density 
matrix e~@( H °~ VH *>\ only the fact that normal-ordered operators are finite when averaged with this density 
matrix. 

Two-loop beta functions. The two-loop result ([50)1 with two-loop divergent parts (fSTj) and one-loop finite 
and divergent parts (169)) and (|73|) is 



- '. A- f - 2 • 2X d (p(w) • In ( A 



10P(w) 



>ln 



- 1 



2P{w) + In 



- 1 



In 



(102) 



This satisfies perturbatively the Callan-Symanzik equation with beta functions 

[3 Xd = -X 2 -X 2 + 0(X 3 ,X 2 X d ) 
f3 x = -2XX d + XX 2 d + A 3 + 0(XX 3 d , A 3 Ad) . 



(103) 



Note that this calculation does not give the third order coefficient of j3\ d . However, we know the universal beta 
functions of the one-channel Kondo model — (g 2 — ^g 3 ) and that of the symmetric two-channel Kondo model 
— {g 2 — g 3 ). Taking V = (which does not affect the beta function), these two cases are obtained respectively 
at A = 0, and at A^ = A upon diagonalization of the matrix of couplings J a , a '- These two facts essentially fix 
the two-loop beta function to be of the form 



p Xd = -(X 2 + X 2 -X 3 -X 2 X d + ...) 
fa = -(2XX d ~aXX 2 -(2-a)X 3 + ...) 

where a is a non-universal number. Our results fix a = 1, which gives the standard beta functions 

Px d = -(A2 + A 2 -A 3 -A 2 A (i + 0(A rf ,A 2 A2,A 4 )) 



fix 



-{2XX d - XX 2 , - X 3 + (9(AA 3 , A 3 A rf )) . 



(104) 



(105) 



It is convenient now to change variables to A± = X d ± A so that (to the same order) 
The RG invariant anisotropy parameter can be expressed as: 



-xl + -x ± (x 2 + + x 2 _) . 



(106) 



(107) 



The scaling limit of the current. We now evaluate the current as function of the voltage and of the 
temperature. For the expression of the current, our analysis will make use solely of the one-loop results ([77)1 . 
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Expression ([77]) is the perturbative expansion in the limit T, V <C A of the model. The couplings Ad, A take 
their "microscopic" values: the actual values for the physical process represented by the interaction terms, and 
the cut-off A is of the scale of the band width. Of course, since we neglected terms vanishing as A — > oo, this is 
not an exact expression for finite A. However, it allows us to have an exact expression in the scaling limit. 

As explained earlier, we expect universality for T, V <C A. It is convenient therefore to consider the limit where 
A is sent to infinity with the couplings becoming cut-off dependent: they are modified so as to keep the physics 
unchanged when the cut-off is increased. Denoting for the moment A = Ao the physical value of the cut-off 
(of the order of the band width) and A = A , A^ = X d the values of the coupling there, the running couplings 
A r (A), A^(A), are governed by the RG equations, 

ra 2 +m 2 -ra 3 -(^ 

2X r X d -X r (X r d f-(X r ) 3 (108) 

with initial conditions fixed by the microscopic values of the couplings: X d (A = Ag) = X d , A r (A = Ao) = A . 

The solution of the RG flow can then be described by RG invariants - quantities that describe the full trajectory. 
Such quantities are C introduced earlier, and a scale Tk to be discussed below. Thus any physical quantity F 
will depend on the cut off and coupling via these invariants F = F(T/Tk, V/Tk, C). One may reformulate this 
scaling procedure as follows. Introduce the scale M , the physical scale on which the system is examined, 

V = M sin(a) , T = M cos(a) 

for some angle a in the V — T plane so that the previous result for the current is written as, 

(J)ss = -7f-V\ 2 

Following the previous considerations, the current can be written in terms of couplings that depend on the 
physical scale M, satisfying the equations with respect to M, 

(AS) 2 -(A r ) 2 + (AT + (A^) 2 AS 

2X r X r d + X r (X r d ) 2 + {X r ) 3 (109) 
with initial conditions X d (M = A) = Ad, X r (Al = A) = A (A being of the order of the band width), as follows: 

{J) ss ^-V{X r f[l + AX r d Q( a ) + ...\ (110) 

where 

Q{a) = P(tan(a)) + ln(csc(a)) . (Ill) 

Again, the solution to the RG flow (1109j) should be described solely as a function of M/Tk and of the RG 
invariant C (|107[) . instead of A and the initial conditions Xd and A. With such a description, we can trivially 
take the scaling limit: T, V, Tk <§C A with fixed ratios T : V : Tk and fixed C, since A does not appear anymore. 
In this limit, the quantum field theory gives exact results, and the system is in its universal regime. In order to 
have unambiguous results, we need to define Tk- Once Tk is defined, one need only solve the RG flow (|109[) 
(a numerical solution is easy to obtain with good precision, for instance), and one obtains the scaling limit in 
its perturbative region. The actual values of Xd, X and of the band width A should be such that the theory is 
near to the scaling limit if we want the system to be meaningfully described by quantum field theory. 

The perturbative region of the scaling limit is Tk -C y/V 2 + T 2 , for any value of C. Hence, in order to define Tk 
in perturbation theory, we need to look at the expansion as M/Tk — > oo of the solution to the RG flow (|109[) . 
This expansion, and the value of Tk, have different forms depending on C . We will consider here two cases: 
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taking C = oo then M/Tk — > oo; and keeping C finite then taking M/Tk — > oo. The first case corresponds to 
Arf = A. In fact, it is true to all orders that if Ad = A, then X r d — X r for all scales M. In this case, the equilibrium 
Hamiltonian H\y=o is in fact a one-channel Kondo model plus a decoupled free massless fermion. The case C 
finite corresponds naturally to A ~ CX d in the scaling limit (as will become clear below). 

In the case where X d = A, the RG equations lead to the following large- (M/Tk) expansion of the running 
couplings: 

_ A , r _ 1 lnln(M/r g ) a In 2 ]n(M/T K ) (4a - 1) In In (M/Tk) / 1 \ 

d ~ ~ 2 \n(M/T K ) + 4 In 2 (M/T K ) + 2 In 2 (M/Tk) + 8 In 3 (M/Tk) + 8 In 3 (M/Tk) + \ In 3 (M/T K )J ' 

(H2) 

The coefficient a can be changed by a change of the scale (the integration constant) Tk', making the replacement 
Tk <— > xTk is equivalent to doing a > a + ln(a;). Note also that a change of scale Tk > xTk corresponds simply 
to a perturbative change of the running coupling constants that keeps the beta functions invariant. Fixing the 
value of a is making a choice of definition for Tk, which is one more condition necessary to totally specify the 
renormalization procedure. Of course, different definitions of Tk reproduce the same scaling limit. For arbitrary 
a, in terms of the coupling A^ = A and of the scale of the band width A, Tk has the form 

T K (a) = A V2\e- a ~ * (i + 0(A)) (A d = A) . 

A standard definition for the Kondo temperature is to make the term in 1/ In 2 (M/Tk) vanish (a = 0), giving: 

T K = A\/2Ae-^(l + 0(A)) (X d = X) . (113) 



In the case where the RG invariant C is finite, the large- (M/Tk) expansion has the form 

= 1 \n\n(M/T K ) a In 2 ]h(M/Tk) (2a - 1) \nln(M/T K ) 

d ln(M/T K ) ln 2 (M/T K ) ln 2 (M/T K ) \n 3 (M/T K ) ln 3 (M/T K ) 
A r 1 21nln(M/TRr) 1 + 2a 31n 2 ln(M/T K ) (6a + l)\nln(M/T K ) 

~C ~ \n 2 (M/T K ) + \n 3 (M/T K ) + \n 3 (M/T K ) + \n(M/T K ) + ln 4 (M/T K ) 



(114) 



where the dots (. . .) mean 0(ln~ 3 (M/TR-)) for X r d , and 0(ln" 4 (M/Tfi:)) for A r . Note that in general, C will 
not appear only as a normalization of A; this is an artifact of the limited perturbative order which we consider. 
Also, note that C must be positive for A r to be positive. Again, a variation of Tk has the effect of changing a: 
making the replacement Tk | — > xTk is equivalent to doing a i— > a + ln(x). For arbitrary a, we can use the first 
or the second equation of (| 1 14[) in order to determine Tk in terms of the couplings A^ , A and of the scale of the 
band width A. This gives two equivalent forms: 

T K (a) =AX d e~ a ~^(l + 0(X d )) = A yAe-o-l-N/? (1 + 0(VX)) (X d ? A) . 

Recall that from p07p . we indeed have that A oc A^ in the scaling limit, so that the second equality above is 
correct. A more standard way of writing the Kondo temperature can be obtained by considering the linear 
combination X d + X r . This gives, again for the same object Tr-(o), 

T K (a) = A(X d + X)e- {a+C) -^(l + 0(X d )) (X d + A) . 
A standard definition is a — — C , giving 

T K - A(A d + A)fT^x(l + 0(A d )) (X d + X) . (115) 



Four comments are now in order. 

First, note that taking C — > oo in (|114p does not give (|112p . This is expected, since the former corresponds to 
taking first M/Tk —* oo then C — > oo in the solution to the RG equations, whereas the latter corresponds to 
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taking first C — > oo then M/Tk — > oo, and these two limits do not commute. In particular, one would have 
obtained still different expressions taking simultaneously C — > oo and M/Tk — > oo, with a prescribed ratio 
between them, for instance (it is a simple matter to evaluate in this case, or in any other case, the expansion as 
M/Tk — * oo of the solutions to the RG equations). The expressions (|114[) and pi2[) should not be regarded as 
describing all solutions to the RG equation; they are rather particular limits of the solutions, as described above, 
and are used here solely for the purpose of defining the Kondo temperature purely using perturbation theory 
(the perturbative regime is, again, the regime M ^> Tk, for any C taken to infinity or not, simultaneously 
or not with M/Tk)- The expansion (|114p offers a way of defining Tk in the regime with C finite, whereas 
the expansion (|112|l offers a way of defining it in the regime with C = oo. Other definitions would have been 
possible, but we only need these two definitions of the Kondo temperature here. 

Second, it is important to recall that the cases C — oo and C < oo do exhaust all possible scaling regimes, 
even though the expressions (|112p and (|114|) do not exhaust all possible behavior of the running couplings as 
M/Tk - > oo- However, there are many ways of reaching any given scaling regime. In particular, the relations 
A = Ad and A ~ CA 2 , (naturally associated to, respectively, the scale definitions (|113|) and (| 1 1 5|) ~) do not exhaust 
all possible ways the same scaling limit can be reached. In order to understand what this means, first recall that 
the full RG-improved perturbation theory reproduces the divergent and finite part of the full bare perturbation 
theory, so that we can talk about the RG trajectory for the bare couplings before taking the scaling limit. 
Then, from this viewpoint, one can take the scaling limit by fixing a trajectory C, and by sending A — * oo 
while keeping A^ and A on the trajectory at scale A. For finite C, this gives A ~ CA 2 , whereas for C = oo, this 
gives A = Xd- But one could also take the scaling limit by changing the value of C (changing the shape of the 
trajectory) while A — > oo, always keeping A^ and A on the trajectory defined by C. For instance, this is what 
happens if one takes A = q\d for some fixed q ^ 1 when sending A^, A — > 0; then one must simultaneously take 
C — > oo in order to keep A and Xd on the RG trajectory. The resulting value of Tk in terms of such couplings is 
different from (| 1 1 3[) and (|115[) . But in this example, when the scaling limit is reached we have C — oo, so that 
the quantum field theory is the same as the one obtained by taking A = A^ and sending them to 0; in particular, 
in the scaling limit, we still have A r — X dl and we can still define a scale Tk using (|1 12[) . Since we are only 
interested in the scaling regimes, we do not need to look at all possible ways a given regime can be reached. 

Third, it is important to note that the leading behavior of the current at large M/Tk is very different if C = oo 
or if C < oo (again, we only look at the two cases mentionned above). In the first case, it is given by 



, on ivi smiLti on v ,. , „. 

~ — ; 2,„L\ = T i 2, /„, . r^, — = A ) • ( 116 ) 



3tt Msin(a) _ 3?r V 

~8 In 2 (M/Tk) ~ ~8~ ln 2 (VV 2 + T 2 /T K ) 

On the other hand, in the second case it is given by 

ottC 2 Msin( Q ) _ 3ttC 2 V 

W)ss 2 \n\M/T K )~ 2 ]n\VW+T^/T K ) ^ ^ ' ' 1 ' 

Note that since the RG invariant C appears as a coefficient, the leading behavior in this case is non-universal; 
this is a property of the regime C finite. 

Finally, the parameter a can be seen as parametrizing a family of choices of infrared cutoffs for our theory. For 
instance, at a = 0, the temperature is the infrared cutoff, whereas at a = tt/2, the voltage solely plays the role 
of an infrared cutoff. It is important to note that the voltage is a good infrared cutoff for the average of the 
current operator to two-loop order. Indeed, the limit T/V — > of our one-loop and two-loop bare perturbative 
results is finite; equivalently, the limit a — > ir/2 of the renormalized perturbative results is finite. We should 
remark however that this need not be the case when other quanties are considered, see Q- 

Universal ratios. From the viewpoint of the interpretation of the measurement of non-equilibrium quantities, 
the Kondo temperature, as defined for instance in (| 1 13|) . is not very convenient, since the temperature T and 
the bias voltage V have a different influence on the Kondo screening cloud, and we would like a definition that 
embodies this difference. Hence, it seems appropriate to define a continuum of scales, Mk(oc), depending on the 
angle a on the V — T plane. Consider first the case Xd = X. A possible definition of Mk(cx) is the requirement 
that the average current J in the steady state does not have a term of the form M/ In (M/Mk (a)) in its 
expansion at large M /Mk{ol)- 

3^ Msmja) f lnln(M /M K (a)) ( \n 2 HM/M K (a)) \\ ( . 

[ ' 8 \n 2 (M/M K (a)) \ ln{M/M K (a)) \ \n 4 (M/M K (a)) J J 
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This gives: 



where Tk is defined in (jTT3j) . 



M K {a) =T K e Q ^ 



(119) 



As particular cases, we can now define the "decoherence" Kondo scale, Mk (tv/2), as the voltage at T = at 
which the Kondo cloud is destroyed by the electrons passing through, and the "thermal" Kondo scale Mk(0). 
We find that their ratio, which is universal and does not receive corrections from the two-loop (or higher-loop) 
contributions to the beta function (nor from higher order contributions to the current), is 
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current 



M K (0) - Z7Te 
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— = 2.96188723. 



(120) 



In the case Xd ^ A (that is, for C < oo), a similar definition of Mk(cx) can be made: requiring that the current 
docs not have a term of the form Mj hi 4 (M / M k (a)) ■ Constructing again the ratio Mk{tt/2)/Mk{0) gives 
exactly the same number: this ratio is indeed universal, independent of both Tr and C . It is important that 
the current possess an infrared-convergent (convergence at large switch-on times) perturbative expansion to 
one loop for these universal ratios to have a meaning. Further, that they have the usual one-loop logarithmic 
accuracy near to the scaling limit is a consequence of the infrared convergence of the two-loop perturbation 
theory. 

The scale Mk{o) defined above is characteristic of the current; other physical quantities would give different 
functions Mr- (a), and different ratios. For instance, the same analysis can be applied on the differential 
conductance 



From the perturbative calculations, we have 



T 



- -1 



Again, in terms of running couplings, we have 



G = — (X r ) 2 
2 v ; 



1 + A\ r d Q{a) 



(121) 



(122) 



where 



Q(a) = P(tan(a)) + ln(csc(a)) + tan(a)P'(tan(a)) 



(123) 



Taking again Xd — A, we can repeat the calculations above and define similarly the scale Mk(ol) associated to 
the conductance. We find 



M K {a) =T K e$W . 
We observe that Q(0) — Q(0) and that Q(n/2) = Q(w/2) — 1. Hence, we have 

Mk{tt/2) 



R, 



conductance 



M K (0) 
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current 
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1.08961742. 



(124) 



(125) 



V. PERSPECTIVES 



Reaching the steady-state. We showed that the bath of free massless fermions suffices to allow the system 
to reach equilibrium at temperature T in the case of zero bias voltage, and to allow it to reach steady state 
when the bias voltage is non-zero. No other relaxation process has to be assumed; the infinite bath of free 
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massless fermions plays the role of a thermal bath and is able to absorb the energy necessary for relaxation 
to occur,as well as for the steady state to exist. This is in close connection with the study of Caldeira and 
Leggett Q: they constructed a model where an infinite number of oscillators provides an explicit dissipation 
in order to study, from first principles, the effect of dissipation on quantum tunneling. It turns out [5( that 
their construction is simply related to models of field theory that are conformal in the bulk (in general, with 
non-conformal boundary conditions), as recall in the Introduction. It would be interesting to see to what extent 
our proof can be generalized to the study of more general models of quantum kinetics with thermal dissipation, 
and to see whether similar arguments can be applied to other impurity models out of equilibrium. A principle 
to follow, as can be extracted from our derivation, is the fact that the orbit of a global symmetry of the model 
should cover the possible values of the boundary degrees of freedom. 

Related to the latter point, another question is: What is the effect of a magnetic field on the steady state in 
the quantum dot? Our proof of factorization in Appendix [B] does not hold anymore when a magnetic field is 
present, since it uses heavily the invariance of the correlation functions under SU(2) transformations. Hence, 
the real-time perturbation series is no longer expected to be convergent as the switch-on time is sent to minus 
infinity. This is simple to understand physically: at small couplings, both the current through the dot and its 
interaction with the leads are weak, that is, both the non-equilibrating effect and the thermalization effect are 
weak. At zero magnetic field, since all states of the isolated dot have the same energy, the thermalization effect 
is more efficient and still stronger. At non-zero magnetic field, however, as the couplings are sent to zero, we 
cannot expect that the dot smoothly reaches its thermal equilibrium energy distribution when no other thermal 
bath is coupled to it. The real-time perturbation theory should describe this situation, but obviously its zeroth 
order cannot give anything else than the thermal equilibrium value of any quantity under study. Hence, in non- 
zero magnetic field, we can expect some strong non-analyticity in the couplings and we must find large-time 
(IR) divergences in the perturbative coefficients. This clearly indicates that the real-time perturbative series 
does not properly describe the approach to the steady state, neither the steady state itself, of the model in 
magnetic field without external thermal bath. 

Questions remain, as raised in the Introduction: Does the quantum field theory still (non-perturbatively) reach 
a steady state, or does it show other behaviors at large times, like oscillations (of the dot magnetization, for 
instance)? If it reaches a steady state, is it a good description of realistic systems, where the dot is coupled 
to an external thermal bath at all times (without exchange of particles), independently of its coupling to the 
leads? 

The former question was partially answered, in [12j | : it was assumed on physical grounds that the non-equilibrium 
Kondo model reaches a steady state, and mainly from this assumption, it was explained how to obtain the zeroth 
order of perturbation theory for the dot magnetization. Indeed, strong non-analyticity is obtained. In a sense, 
one should start with a density matrix that already contains a non-thermal distribution of the impurity spin 
states; this density matrix can be obtained by requiring that the perturbative series be convergent at large times. 
As explained clearly there, this is equivalent to solving a quantum Boltzmann equation in order to determine 
the non-thermal dot occupation numbers. This should answer partly, in some sense still perturbatively, the 
question of describing the steady state of the quantum field theory (the non-equilibrium Kondo model with 
magnetic field) - although the results really start with the assumption of a steady state, and do not establish 
its existence. However, this does not address the question as to whether the leads correctly play the role of 
thermal baths, or whether a coupling to an external thermal bath would have important effects. 

In relation to the latter question, our argument suggests that, at least at small couplings (or at temperatures 
much greater than the Kondo temperature), the model does not describe the true steady state of the non- 
equilibrium Kondo dot in contact with a thermal environment. Indeed, from a physical interpretation of the 
perturbative series, the leads do not provide a strong enough thermalisation to absorb the energy necessary 
for the steady state to occur, and it is possible that it cannot be trusted to sustain the correct steady state. 
That it does not provide the thermalisation for the steady state to occur is certainly in agreement with [Tig : 
there it was one of the main points that the large-time divergences are due to the absence of a proper thermal 
bath, and that one needs to put "by hand" a thermal bath connected to the dot. This was done, essentially, by 
putting a small imaginary part on the evolution time, which was then set to zero before taking small couplings 
A, Xd (that is, it was obtained a steady state where the coupling to the thermal bath is much smaller than the 
couplings A, Ad). The most delicate question, however, concerns the fact that the steady state itself may be 
affected non-trivially by a thermal bath. It is probable (but this should be verified) that other, more realistic, 
representations of a thermal bath give the same results as [l2| in the limit of small coupling with the thermal 
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bath. However, as the couplings A, A^ are sent to zero, due to the thermal effect of the environment, the average 
dot magnetization, for instance, should smoothly reach its thermal equilibrium value, and this may well be a 
universal cross-over behavior (since it occurs at small couplings, hence near to a second order phase transition). 
In order to assess this, it would be important to have a more adequate description of a thermal bath, and the 
theory of Caldeira and Leggett surely provides the most promising avenue. 

The steady-state density matrix. Our proof that the real-time perturbation theory describes a steady 
state was immediately adapted to the proof that the steady state can be described by a density matrix (140)) , 
almost as in an equilibrium state. As mentioned, a steady-state density matrix was also introduced in [6| for 
generic models under the assumption that a relaxation time was present. Our derivation is slightly different, 
and does not make further assumptions, since we showed that relaxation does occur in the non-equilibrium 
Kondo model. The main characteristics of the steady-state density matrix, as opposed to equilibrium density 
matrices, is that it is defined by coupling a non-local conserved charge Y to the voltage V, instead of coupling 
a local conserved charge to an appropriate chemical potential. This non-locality is related to the fact that the 
operator Y describes the build-up of the steady state, or, in some sense, the asymptotic state that characterizes 
the steady state. The properties of this asymptotic states and other consequences of this description are still 
to be explored. 

Renormalized real-time perturbation theory. We developed the two-loop renormalized real-time pertur- 
bation theory. We gave an argument for the validity of the Callan-Symanzik equation in the steady state (with 
the same beta functions as the equilibrium ones), and verified this to one loop for the current operator inside 
any correlation function. It is tempting to relate the validity of the Callan-Symanzik equation to the fact that 
the operator Y appearing in the steady-state density matrix is a conserved charge. In particular, the fact that 
the voltage does not flow is obvious from such consideration, as it is coupled to a conserved charge. However, 
since Y is non-local, it is hard to make this connection more precise. 

The quantum field theory gives physical quantities in the scaling limit (the universal region) V, T, Tk <C D, 
where Tk is the Kondo temperature, which is a non-universal quantity related to the microscopic values of the 
couplings. Our renormalized perturbative results give the current in the region Tk <C y/V 2 + T 2 <C D. This 
includes the part of the universal region where the system is strongly out of equilibrium. In this perturbative 
region, we defined a continuous family of Kondo scales M(a) depending on the ratio V/T = tan(a). By 
comparing the scale at a = tt/2 (T/V — ► 0) with that at a — (V/T — + 0), we obtained a universal measure of 
the effect of the voltage on the Kondo cloud, as compared to the effect of the temperature. We noted that such 
a universal measure was correct since to two-loop order, the voltage plays the role of a good infrared cutoff for 
the current, so that no divergencies appear to that order as T/Tk — > (with fixed V/Tk)- 

In connection to the latter point, it has sometimes been suggested in the literature that at T /Tk = 0, the system 
should be in a "strong coupling regime" (see for instance [18J), and as such, the perturbation theory should 
not be valid and should show infrared divergencies (divergencies as T/Tk — * 0) in higher-loop calculations. In 
particular, it is clear that this occurs in the calculation of any thermodynamical quantities, which can indeed be 
deemed "in a strong coupling regime" (at the IR fixed point) at zero temperature. In [8, 9], for instance, one sees 
logarithmic divergences as T/Tk - > for any fixed V/Tk at the one-loop order of perturbation theory for the 
spin susceptibility, pointing to the fact that for describing the limit T/Tk — > of that thermodynamical quantity, 
one needs to know about the IR fixed point [231] (the perturbation theory only describes the theory around its 
UV fixed point). Our point, though, is that this might not be so for all quantities. Our two-loop results suggest 
that in a sense, the current is really a dynamical quantity, ruled by the scale V. For comparison, this is much 
like a correlation function of two local fields is ruled by the distance between the points in equilibrium quantum 
field theory, even at zero temperature (the short distance behavior is described, for instance, by the UV fixed 
point of the theory, up to a normalization if the fields are not conserved currents, and up to the one-point 
functions of the operators appearing in the operator product expansion of the fields). Then, the region V 3> Tk 
should really be, for the current, a weak-coupling, UV situation, even at zero temperature; this is at least what 
we see at two-loop order. Notwithstanding the fact that infrared divergencies may signal that the steady state 
is not reached, as discussed above, possible divergencies in the perturbative expansion of the current at higher 
orders may correspond to simple power-like non-analyticity in the couplings at zero temperature, with smaller 
contributions than those of the two-loop results (as in the usual situation of correlation functions in equilibrium 
quantum field theory). 

Finally, it would be very interesting to fully verify the validity of the Callan-Symanzik equation with a magnetic 



32 



field. It is easy to check for instance that the one-loop, finite-magnetic-field results of @, S El satisfy the 
one-loop Callan-Symanzik equations in the universal regime. Although we did not cover the case with a finite 
magnetic field (and see the discussion above for the subtleties involved) , we expect that our general arguments 
for the validity of the Callan-Symanzik equations to all order still hold since in the universal regime, the magnetic 
field is a low-energy scale as compared to the band width. It may be useful to note that the usual perturbative 
renormalization was modified in [H, [^] in order to correctly incorporate the structure of logarithmic divergencies 
in the region V > D (where D is the bandwidth) of the one-loop calculation of the current and of the impurity 
magnetization at finite magnetic field; more precisely, energy-dependent coupling constants were introduced. 
We want to stress that this region may be non-universal (that is, results may depend on the precise structure 
of the band), as is the region T > D or Tk > D. Then, naturally, it cannot be universally described by a finite 
number of coupling constants (or, more precisely, by a finite number of RG invariants). In order to recover 
"scaling" , one needs to use exact RG or similar methods, and the usual Callan-Symanzik equation does not 
hold; but this is not in disagreement with our results, which only deal with the universal regime V, T, Tk <C D. 
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APPENDIX A: CONNECTED CORRELATION FUNCTIONS 



Consider a theory with density matrix p, and denote, for any operator O, 

Tr [pO] 



((O)) 



Tr[p] 



Consider the following average: 



((^xp(V t >lMu))) 

where Hi(t) can be any operator depending on the (time) parameter t, and O is also any operator. The path- 
ordered exponential above is understood as an expansion in time-ordered integrals of multilinears of Hj(t)'s. 
The connected correlation functions, where the Hj(t)'s are "connected" to O, can then naturally be defined by 
saying that 



oo „( 
„=0 Jta 



• 

dh I dt 2 --- I dt n ((H I (t 1 )H I (t 2 )---H I (t n )0)) conncctcd . (2) 



Equivalently, they can be defined recursively by 

((Hj(t 1 )---H I (t n )0)) 

= ((H I (t 1 )---H I (t n )0)) 

connected (3) 
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This recursive definition does not involve integrations over the time parameters tj's, so that it is in fact a general 
definition for correlation functions where a set of arbitrary operators i?/(tj)'s parametrized by i = 1, . . . , n are 
connected to O. 

The main property of connected correlation functions is the following. Consider a fixed set of U's for i = 1, . . . , n 
and a fixed subset t ai , . . . , t am for m < n. Consider correlation functions of i2j (ij)'s with and without insertion 
of O. One can observe that if the operators Hj(t ai Ys factorize out of all such correlation functions (for instance: 

{{Hi(h) ■ ■ ■ Hj{t ai ) ■ ■ ■ H T (t am ) ■ ■ ■ H!(t n ) O)) 

= {{H I (t 1 )---H^ 1 )---H^(tZ)---H I (t n )0}) ((H I (t ai )---H I (t a J)) ) 

then all connected correlation functions involving at least one of the Hi(t ai Ys are zero. This is easy to show 
from (J21) by induction on the number of operators Hi(tiYs inside connected correlation functions. 

Connected correlation functions also occur in more general situations: 

'v exp U fi! dt H { j 1] (t)) V exp (i fif dt Hf ] {t)Y--V exp (i fif dt H { j N) (t)j O 



V exp (i f*} dt (*)) V exp (i fi* dt H$ ] (t)) • • • V exp (i fig dt h\ N) (t) 
= //vcxp (iJdtHP(t)\ Pcxp (i JJltHf\t) \ ■'•VexpfiJ dtH$ N) {t) \ O 




connected 



On the right-hand side, the operators H\ (ti)'s, Hi (ti)'s, H\ (t,)'s are all connected to O. 



APPENDIX B: PROOF OF FACTORIZATION 



Consider a product of operators of the type 

Ji(x + xi) ■ S J 2 (x + x 2 ) ■ S ■ ■ ■ J n (x + x n ) ■ S 

where Ji,2,... can be J<j, J x or J y . Recall the notation (|19| . There, the trace is performed over the Hilbert 
space of the conformal field theory where the currents Ji,2,... act and over the two-dimensional impurity space 
associated to S. Insert an operator which is composed of products of local operators acting on the CFT Hilbert 
space tensored with an arbitrary operator on the impurity space. We will denote it by 

a 

where a = 0,1,2,3 and Sq is the identity 1 on the impurity space. The operators in O a can be at any fixed 
time (with respect to the theory Hq) and position, and can also be integrals of such operators over finite time 
intervals. We the consider the SU(2) invariant quantity: 

«Ji(z + Xl )-S J 2 (x + x 2 )-S--- J n (x + x n )-S O a S a )) 

= ((Jl 1 (x + x 1 )J?(x + x 2 )---.P n "(x + x n )O a )) ((S ll S l2 ---S ln S a ))o . (1) 

The first factor on the right-hand side factorizes as \x\ — > oo because of the locality of the operators and because 
the correlation function is evaluated in a unitary quantum field theory: 



((j^(x + Xl )j^(x + x 2 )---r r :(x + x n )o a ))o ~ ((j?(x 1 )j?(x 2 )---r r r(x n )))o((o a ))o (2) 

This field theoretic factorization induces a corresponding factorization in the spin space, as we proceed to show. 
The expression ((^(xi) J^{x 2 ) ■ ■ ■ J^"(x n )))o is a tensor in the product space of n copies of the fundamental 
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representation of rhe rotation group 0(3) and the only way to form a rotational invariant with it is to multiply 
it by a corresponding product of matrices S^S^ ■ ■ ■ Si n , i.e. (( J] 1 (x%) J % 2 X (x 2 ) • ■ • J„" (x n )))o Si ± Si 2 ■ ■ ■ Si n oc 1 . 
Hence, the product S{ 2 ■ ■ ■ Si n also factorizes, and we have 

«J^i) J£(x 2 ) ■ ■ ■ J£(x n ))) <(CT)}o ((S h S i2 ■ ■ -S in )) ((S a )) Q (3) 
The factorization property then follows: 

{{Jx{x + x x )-S J 2 (x + x 2 )-S---J n {x + x n ) -SO a S a )) 

N ~°° ((Mx + Xl ) ■ S J 2 (x + x 2 ) ■ S ■ ■ ■ I n (x + x n ) ■ S)) ((O a S a )) Q . (4) 

The operator Hi is a linear combinations of operators of the type J ■ S. The operators J are right- moving fields. 
Hence when Hi is evolved in time with Hq for a time t, it becomes a linear combinations of J • S"s at position 
—t. This shows the factorization of juxtaposed Hj (t)'s as \t\ — > 00 when they are evaluated inside a trace with 
operators at fixed time and position on the right. The same proof applies if such operator insertion is put on 
the left of juxtaposed ij[°^(t)'s. This completes the proof. 



APPENDIX C: CALLAN-SYMANZIK EQUATION FOR MATRIX ELEMENTS 



In this appendix, we will justify equation (|87[) . Consider the finite-temperature average, in the theory described 
by -ff |y=o, of the regularized current operator J7a = A(J x )a(0) • S with insertions of creation and anihilation 
operators of the Hamiltonian: A\{p) and Ai(p). More precisely these are creation and annihilation operators 
for eigenstates of H \v=o with energy p, corresponding to the massless particles "naturally" associated to the 
local operator Ji(x), where Jj is any of the ten operators Jd, J x , Jy, Jz- Thus consider, 

Tr ( e -WW=o Ail (pO • • • A im ( Pm )J A A[ ( Pm+1 ) ■ ■ ■ A\ ( Pm+n )) 

\ m+l m+« )_ . , 

Tr ( e -^ ff lv=o) ■ 1 ' 

Since the creation and annihilation operators are eigenoperators of the Hamiltonian, the quantity ([1]) satisfies 
the Callan-Symanzik equation (|83|) (that is, in ([83| we can replace the steady-state average of the current by 
this quantity). 

But we showed that the interacting density matrix e _/3i/ ' v '=° can be obtained from the free one e~^ H ° by 
evolving it for an infinite time (| 1 5 j) . In much the same way, the interacting creation and annihilation operators 
can be heuristically written in terms of the mode operators associated to the current algebra operators: 

A\(p) = S(Q,-oo)4(p)«S'(-oo,0) , Ai(p) = S(0,-oo)Ji(p)S(-oo,0) 

(where S(t\,t 2 ) is defined in Q). Here, Ai is any of the 10 operators Ad, A x , A y or A z . Indeed, the operators 
A\(p) and written in this way heuristically satisfy the canonical commutation relations amongst them, 

and the appropriate commutation relations with Hamiltonian H\y=a, written as >S'(0, — oo)HqS(— 00, 0) (they 
are eigneoperators of H\v=o)- Hence, the quantity |1| can be written 

Tr (e-^J^Cpi) • • • Ji m (p m )S(-oo,0)J A S(0,-oo)jj m+1 {p m+1 ) ■ ■ ■ 4 m+ „(Pm+n)) 

Taking the zero-temperature limit /3 — * 00, this quantity reproduces all matrix elements of the current (see 
ean ([5S|) ) Ja = S(—oo, 0)JaS(0, —00) on the Hilbert space of the free theory Ho, which proves ((57)) . 

In fact, to be more precise, when going more from to Jl}, one needs to modify slightly ©. First, the mode 
operators jj{p), Ji{p) should be replaced by appropriate wave packets jj{p), Ji(p), obtained by integrating 
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the associated currents Ji(x) (times the oscillating exponential e~ lpx ) over space with a kernel vanishing (ex- 
ponentially, say) as \x\ 3> I for some wave-packet extent I. Then, the evolution operator 5(— oo, 0) should be 
replaced by S(io,0) and 5(0, oo) by 5(0, t' ), and these two evolution operators should be put into the trace in 
the denominator as well, in the order in which they appear in the numerator. By using arguments involving the 
vanishing of connected correlation functions, it should be observed that the limits |to| 3> 1,(3 and \t' \ 1,(3 exist 
independently. Bringing both evolution operators around the density matrix e~^ H °, one should then use the 
steps involved in proving (|15p . The result, up to vanishing contributions as to,t' 3> 1,(3, is (|TJ) , with operators 
A\(p) and Ai(p) replaced respectively by 

A\{p) =S(0,t' +i/3)jj(p)S(t ,0) , Aiip) = S(0,t' + H3)Ji{p)S(t ,0) . 

These are the operators that should correspond to asymptotic states, in the limit |io|, |£q| ^ I ^ (3 with 
to = t' < 0. In this limit, one indeed recovers @. 



APPENDIX D: ASYMPTOTIC BEHAVIOR OF THE FUNCTION P(w) 



In this appendix we evaluate the asymptotic behavior of the function P(w) ([74]) . The large w behavior is easy to 
obtain; the first term inside the parenthesis (on the right-hand side of the second equation in (|74p) disappears, 
and the second term — (p — l)/(e" , ' p_1 ' ) — 1) becomes (p — l)8(p < 1): it is non-zero only for p < 1. Taking 

2 

these contributions into account as well as the last term e~ p , the result is: 

P(oo) = l-|. (1) 

In fact, we can also obtain the next term in the large w expansion. The next contributions can be written in 
the following way, by expanding the integrand: 

/•OO J 

/ — y f (P + l)e~ nw[p+1) + 0(p < l)(p - l)e nw ( p ~V - 0(p > i)(p - iy^-p) 

Interchanging integration and summation (this is valid because the integration variable p is always kept in 
a region where the expansion of the integrand is convergent), the integral of every term of the sum can be 
expressed in terms of the exponential integrals. Every term can then be expanded at large w, giving 

^ / 2 12 240 

' V n 2 w 2 n 4 ui 4 n 6 u) 6 

n— 1 v 

Now every term can be re-summed, and we get the large- u> asymptotic expansion (|75p . 
The small w behavior is more subtle. The integral can be divided into two parts: 

The first integral can be evaluated by expanding the integrand in small w. The integrand goes as — l+e~ p " +0(w) 
at small w, so the first integral is convergent at w — > 0, and in fact gives an expansion in Taylor series in w. In 
order to obtain the divergent part in w, we need only consider the second integral, and we can forget about the 
term e~ p . Make the transformation of variable p — > p/w: 

dp ( p/w + 1 p/w — 1 



p \ eP+ w - 1 eP~ w - 1 

The integrand can then be expanded in w: this gives a Taylor series in w 2 starting with power 0. Each term of 
this Taylor series gives a convergent integral: the asymptotic behavior of each term at p — > oo is exponentially 
decreasing. Moreover, each term, except for the very first one, has a behavior like p° as p — > 0, so that at each 
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order in w 2 , except at the zeroth order, we can evaluate the integral. We thus obtain a Taylor series in w 2 . At 
the zeroth order, we have 

dp ( 2 ( P -l)eP + l 



(eP - l) 2 



We can write it as 



dp ( _ 2 (p~l)e p + 1 



p V ( eP 



P T ±l + e(p< i)) + £° ^(-e(p < i)) 



The first integral is convergent as w — > 0: it has again a Taylor expansion in w. The second integral is easy to 
evaluate, and gives ln(ui). Hence, 

P(w) ~ ln(w) as w -> . (3) 
In fact, we can gather the previous missing parts to get the constant term: 



This gives 



so that we get ([To]) . 



^ +7- ln(27r) = -0.374434476... (5) 



APPENDIX E: INTEGRALS FOR THE TWO-LOOP CALCULATIONS 



The integrals, as they enter in ([50)1 . are 



/ <ir R{r) ( -inS(p + q + r) + P ) R(-p - q - r) x 

7 V p + q + rj 

( - 2 ff (p + q)(-f(p + q + r)+ f(r)) + g(q)(-f(p + q + r) + /(p + r) - /(g + r) + /(r)) 



h = -^f- J d P R(p) f-i n 5(p)+P^\ J dqR(q) (-inS(p + q) + P^ 

drR(r) \-inS(p+q + r) +P ) R(-p-q-r) x 

V ' p + q + rj 



( - 2g(r)(f(p + q)-f( P + q + r)) + ff(r)(/ («) - /(g + r)) + g(p + r)(/(g) - /(p + g + r)) 

-^p | dpi?(p) 2 (p^) | dgi^g) 2 LiTT^ + gJ+P-i-^ (-i*8(q) + P^ p/(g) 
-^p 1 dpR{pf{-m5{p)) y dgi^g) 2 N^(g)+P^Jp/(?) 



^ 1 dpfl(p) f-«rJ(p)+pM y dqR(q) (-^(p + gJ+P-l^ 

y <ir i?(r) ^-ztt(5(p + g + r) + P ^ - * - - j i?(-p - g - r) X 

(2g(r)(f(p + q)-f( P + q + r)) + g(q)(-f(p + q + r) + f(p + r)) + 4) 
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where 



~ J dpR(p) 2 (iw6'(p)+P± 



dqR{q) 2 -ill8(p + q)+¥- 



1 



p + q 
1 



g.g(g) 



I 7 = ~- j dpR{p) ^-iTT^O) + J j dqR(q) r-i^(p+g) + P_l_ 



y dr i?(r) ^-ztt5(p + q + r) + - * - - j i?(-p - g - r) X 

(2 fl (p + g)(-/(p + q + r) + f(r)) + g(p + r)(f(q) -f(p + q + r)) - <?(<? + r)(f(q) + f(r)) + 1 



dgi?(g) 2 -iirJ(p + g) + P- 



-iir5(q) 



P- ) p/(«) 
9/ 



■^p | dpJ2(p) 3 (-i7r*(p)) | dqR(q) 2 (in8'(q)+P^) pf(q) 



R(p) = R A (Vp) = e 



= e 2A^ 



ff(p) 



1 



1 _ e -wp 

Recall that w = (3V. The function f(p) is as in (fTTj) , and is related to g(p) by 

/(p) = (p + l)$(p + 1) - (p - l)ff(p - 1) . 



(1) 



(2) 



(3) 



APPENDIX F: FINITE CONTRIBUTIONS OF THE TWO-LOOP RESULTS 



They are given by 

[h + h + ^6 J finite 



( . t1 . j 00 dq r*dr {gc(q + r){m + /(r)) _ Uq _ r){m /(r))) 



o V Jo r 
3nV f°° dq f q dr 

2 Jo Q Jo r 
3nV [°° dq f q dr 

(l q Jo r 
dqdr 



2 

2 ./o 1 r 
f°° dqdr 

o 9 r 



9c(q - r)f c (q - r) - g c (q + r)f c (q + r) 



M<7 + r ) h c {q - r)e" 

q — r 



(g c (r) + e- r )(2/(g) - f(q + r) - /(g - r)) 
(l-e- r2 )(f c ( q + r) + f c (q-r)) 



UttvJ^ Q + In r)f c (r) 



3nV 
3ttV 



dq 

2~Jo T 

3tt 3 V ( 3 



/c(g) 



+^(-GH-i 

2 \ w 



+37rV ( u(l -i)+u(l + i)-u 

n 3 V 
24~ 



1 - i 



l + i 



(1) 
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and 



[I: 



3ttV 



7j finite — 



Here, we have 



dq 



dr r 



+ 



2 Jo 9 io r 
3ttV r°° dq f q dr 

2 Jo 9 io r 
3irV f°° dq f q dr 

2 Jo ?i r 

- Q nv roar* 

Jo q Jo r 



9c(q - r)(f(q) - f(r)) - g c (q + r)(f(q) + f{r)) 

9c(q - r)f c (q - r) - g c (q + r)f c (q + r) 
Q 



h c (q + r) 



-h c (q - r)e" 



q — r 

(l-e- r2 )(l-e-« 2 )f c (q + r) 



+3ttV 



r°° dr n 

Jo 



+ lnr)/ c (r) 



3irV f°° dr n , \ , / x 



3ttF 



+6irV 



dq 



fc(q) 



u(i — 1) + u(—i — 1) + u 



i - 1 



+ u 



-i - 1 



37i-yin2 5tt 3 V 
+ ^ + —r- 



g(p) 
/(p) 

5c 0) 

hip) 

h c (p) 
G(w) 

u(x) 



1 



1 — e~' u ' p 



(p + l) 5 (p + 1) - (p - l)fl(p - 1) 
50) - sign(p) 
.f(p)-2sign(p)(l- e -P 2 ) 

/ c (p) + (l-e-f 2 ) 3c (p) 
sinh(ui) — ^ 

sinh 2 (f ) 
— dilog(l — x) 



( 2 ) 



(3) 



where to = V/T. 
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